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Abstract 

We study on the lattice the 3d SU(2)+Higgs model, which is an effective theory of a 
large class of 4d high temperature gauge theories. Using the exact constant physics 
curve, continuum (V^ — 00, a — 0) results for the properties of the phase transition 
(critical temperature, latent heat, interface tension) are given. The 3-loop correction to 
the effective potential of the scalar field is determined. The masses of scalar and vector 
excitations are determined and found to be larger in the symmetric than in the broken 
phase. The vector mass is considerably larger than the scalar one, which suggests a 
further simplification to a scalar effective theory at large mn- The use of consistent 
1-loop relations between 3d parameters and 4d physics permits one to convert the 3d 
simulation results to quantitatively accurate numbers for different physical theories, 
such as the Standard Model - excluding possible nonperturbative effects of the U(l) 
subgroup - for Higgs masses up to about 70 GeV. The applications of our results to 
cosmology are discussed. 
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1 Introduction 



The 3d SU(2)xU(l)+Higgs model is a universal theory for the description of the 
electroweak phase transition in the standard electroweak theory and many extensions 
thereof, including the MSSM [1-5] (for a motivation of the study of the electroweak 
phase transition, see [1-4]). In the present paper, we study with lattice simulations 
the dominant SU(2)+Higgs part of the theory, defined by the Lagrangian 

L = -^F^^F^^ + (A0)^(A0) + m20t0 + A3(0V)'. (1.1) 

The procedure of dimensional reduction [0], [6-13] allows one to compute perturba- 
tively the relationship between the temperature T and the physical parameters of the 
underlying 4d electroweak theory or its extensions, and the parameters of the 3d the- 
ory. Concrete formulae for the SU(2)+Higgs model and MSM can be found in [Q (see 
also below). 

The aim of the present paper is to study the 3d SU(2)+Higgs model, especially its 
phase diagram, on the lattice. We confine ourselves to a small ratio X^/g^ < 1/8, which 
in 4d terms corresponds to the case of small Higgs masses, rriH < ttlw ^ 80 GeV. This 
case seems to be the most interesting one for cosmological applications, because the 
phase transition at least in this region is of first order. 

First Monte Carlo results on the electroweak phase transition using a 3d effective 
theory have already been given in [0, |^, see also [0, |15|. In comparison with [|I|, ^ 



we considerably extend numerical calculations. This makes it possible to determine 
for the first time or more accurately than previously a number of finite T quantities 
such as latent heat, correlation lengths, interface tension, magnitude of the higher 
order perturbative terms, etc. An essential ingredient in the increased accuracy is that 
the continuum-lattice mapping formulae, which are exact in 3d, are now known |§, 



T6| . Thus the continuum limit V — > oo,a can be carried out under controlled 
conditions. Preliminary results of the simulations described here were published in [jl^ . 

Lattice Monte Carlo studies of the 4d SU(2)+Higgs model have been reported in [18- 
23]. Whenever comparison can be made, the results are in agreement within errorbars. 
However, the 3d approach used in this paper gives much smaller errors than the 4d one. 
Monte Carlo simulations of the 0(4) pure scalar theory in 3d with a non-analytic cubic 
term have been performed in [|T^]. The spectrum of excitations in this theory, however, 
is very different from that of the SU(2) model, so that comparison is not possible. 

The paper is organized as follows. In Section |^ we present the basic relations for the 
3d SU(2)-|-Higgs model in continuum and on lattice. In Section ^ we describe the Monte 
Carlo update algorithm. In Section ^ we study the properties of the broken phase and 
the convergence of perturbation theory there. Section ^ is a lattice investigation of the 
phase transition for different values of the scalar self-coupling. In Section ^ we measure 
the different correlation lengths, and in Section |^ we study the metastability region and 
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the properties of the symmetric phase. Section | contains a short account of simulations 
with a larger Higgs mass, and in Section ^ we discuss the problem of studying large 
Higgs masses more generally. In Section |10| we summarize the information acquired on 
lattice about the properties of the 3d SU(2)+Higgs theory, and compare to perturbation 
theory and to some non-perturbative approaches. In Section |ll] we relate the 3d lattice 
results to 4d continuum physics. Section ^ is a discussion of some implications of our 
results to cosmology. The conclusions and proposals for future work are in Section p^ . 

Readers not interested in details of lattice simulations could check tables 5-7, and 
go directly to Section [1^, which contains a summary of the non-perturbative results. 

In this paper we use results from and . All references to specific formulae from 
these papers are indicated by I and II followed by the number of the corresponding 
expression. 



2 3d theory in continuum and on lattice 

To make the paper self-contained, we summarize here the essential properties of the 
3d theory in continuum and on lattice. A more detailed discussion can be found 
in §1,0. 

The two couplings and A3 of the theory in eq. ( [I.l| ) have the dimensionality of mass. 
They do not possess ultraviolet renormalization and, say, in the MS scheme are scale 
{n) independent. The mass squared of the scalar field has a linear and logarithmic 
divergence on the 1- and 2-loop levels, respectively. In the MS scheme the relation 
between the scalar mass and the renormalization-group invariant parameter A^ is 



mlin) = ^,log^, (2.1) 
where, in the SU(2)+Higgs theory. 



f2m = ^gt + - 12A2. (2.2) 

Since all the three parameters of the 3d theory are dimensionful, the theory is 
uniquely fixed by giving three parameters, the gauge coupling of dimension mass and 
two dimensionless ratios: 



gl ,.!!i(||).J,|^lo,%. (2.3) 



Renormalization introduces an intermediate mass scale on which the physics does not 
depend. 
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The 3d theory of eq. (|1.1|) is described on a lattice with the lattice constant a by the 
action 

1 



S = /3GEE(l-^TrP,,) + 

X i<j ^ 

- /3/^EE^Tr$t(x)C/,(x)<l>(x + ^) + 

X i ^ 

+ E ^Tr $t(x)$(x) + (5n E[^Tr $t(x)$(x) - 1] 



(2.4) 



The three dimensionless parameters (3giI3hiI3r of eq. ( |2.4| ) are in the continuum limit 
a — *• related to the three dimensionless parameters gla, x, y by the following equations: 
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2^^-^^'i + C)+5.0 + 5.2a; 



(2.5) 
(2.6) 

(2.7) 



Eq. (|2.7] ) depends on several constants arising from lattice perturbation theory; S = 
3.17591, C = 0.09 and the two numbers 5.0 and 5.2, specific for SU(2)+Higgs the- 
ory, and computed in |T^. This reference also gives the analogous numbers for some 
other relevant 3d theories. Note that the logarithmic 2-loop term on the second line in 
eq. ( |2.7| ) is absolutely necessary with the accuracy which we have in our lattice simu- 
lations. For instance, changing the number 5.0 in the 2-loop part by 0.05 changes the 
critical temperature by an amount equal to the statistical uncertainty in one of our 
lattice simulations (138.38±0.05 GeV). 

When the lattice constant a is varied, eqs. (|2.5|-|2.7|) define for the fixed parameters gl, 



X, y of eq. ( p.3|) a curve, the constant physics curve (CPC), in the space of I3gi I3hi Pr- 
All the curves end in the point (cxo, 1/3, 0) for a ^ 0. 

The above discussion was entirely confined to the 3d theory. As explained in [Q, a 
single 3d theory is the effective theory of a large class of finite T field theories and we 
shall later give several quantitative examples of this. For each 4d theory one separately 
has to establish the transformation from the physical parameters of the 4d theory to 
gl,x,y. For the Standard Model these are given in 

Since the use of the 3d parameters g"^, x, y is rather unilluminating, we shall in the 
presentation of lattice results replace them by a "Higgs mass" 171*^ and "temperature" 
T* using the following equations: 



9l 



0.44015r* 



(2.8) 
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X = ^ = -0.00550 + 0.12622/i2, (2.9) 

93 



9l 



(2.10) 



* \2 

rri i ^ 



0.39818 + 0.15545/^2 - 0.00190/i^ - 2.58088 , 



where 



/i^ . (2.11) 

80.6GeV ^ ^ 

These equations follow from the tree-level formulae 

1 im* )2 

g = 2/3, A = -/ ^ , = 80.6 GeV, = W^^T = 0.60858T, (2.12) 




from a formula relating mKfi) to 4d quantities (1.66), and from a subsequent integration 
over the Ao-field 0. Each quoted set of (m^, T*) can with eqs. ( |2.8| - ^.1U| ) be converted 
to a set of gl,x,y, which for a given theory discussed in Sec. |lT], in turn be 

converted to a precise set of values for the zero-temperature pole Higgs mass and the 
physical temperature. For the SU(2)+Higgs theory without fermions the difference 
between the simplified parameters m|^, T* and the physical parameters is relatively 
small. 

In studying the phase structure of the 3d theory we search for a critical curve in 
the (x, ?/)-plane. On the tree level this curve is the line y = 0: for y > the theory 
is in a symmetric, for ?/ < in a broken phase. In 1- or 2-loop perturbation theory 
the line splits in three: a critical curve y = ydx) and upper and lower ends y = y±{x) 
of metastability branches. On the critical curve the system can exist in two different 
phases with the same vacuum energy evac but different values of various gauge invariant 
condensates like The broken phase exists for y < y+(x), the symmetric phase for 

y > y-{x). One knows that perturbation theory can never conclusively determine the 
curve y = ydx) or the jumps of various gauge invariant condensates like across 
the curve, and the main purpose of the lattice Monte Carlo study is to do this. For 
instance, one is interested in knowing whether the curve y = ydx) continues to large 
(> 1/8) values of x or whether it terminates. Our present simulations do not provide 
an answer to the last question. 

When performing simulations with the action (|2.4|) , the procedure is somewhat dif- 
ferent depending on whether one is performing simulations at some fixed {x, y) or 
whether one is searching for the critical curve y = ydx). 

Assume first that one is studying the system, e.g., determining correlation lengths 
or the value of (0^); at some fixed {x,y). The procedure then is as follows. Choose 
some (3g which then gives the lattice constant a. The value of on an A^^ lattice 
should satisfy the constraints that the smallest correlation length is larger than a and 
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the largest correlation length is smaller than Na. The former requirement gives a lower 
limit for f^c and the latter a lower limit for A^. Since (x, y) are fixed, eq. ( |2.7| ) gives the 
value of Ph and eq. ( p.6|) the value of Pr. Simulate the system with these values for 
larger and larger N and perform an extrapolation to — oo. Choose then larger and 
larger values of Pg (smaller and smaller a) doing always the same. The set of different 
extrapolations to — * oo can then be extrapolated to jSc oo (a ^ 0), which is the 
final continuum limit. 

When searching for the critical curve ydx), one can in the beginning only fix the 
value of X. Then one again first chooses Pa subject to the constraints given above, 
simulates the system for various Ph with Pn = xPjj/Pc and finds the value of Ph,c- 
Extrapolating at fixed Pg to N = oo gives a value Ph,c{N = oo) which using eq. ( p^) 
can be converted to a value of y{x) at this a. The extrapolation to a ^ is carried 
out as before and gives the final y = ydx). 



3 The Monte Carlo update algorithm 

The lattice spacing a and the linear lattice size Na are constrained by the length 
scales set by the W and Higgs masses: a <^ l/m;4/(T) < l/niHiT) <^ Na (assuming 
mH{T) < my/{T)). Even though this requirement is much milder than the 4d one 
(a <^ 1/^) 0]) many cases it still mandates quite large lattice sizes (our largest 
volume is 50^ x 200). Therefore, it is important that the update algorithm be as 
efficient as possible. 

The gauge field update is not qualitatively different from the standard SU(2) pure 
gauge update, in spite of the hopping term Tr $'''(x)[/j(x)$(x + i) in the action. To 
update the gauge links we use the conventional reflection over relaxation and Kennedy- 
Pendleton heat bath |Q methods. All the gauge field modes are much 'faster' than 



the Higgs modes, i.e., they have much shorter autocorrelation times. 

Due to the flatness of the Higgs potential, the 'slow' modes of the system are asso- 
ciated with the radial sector of the Higgs field ^ = RV , R>Q,V & SU(2). In what 
follows we summarize the more non-standard methods we use to increase the efficiency 
of the update program. 



Global radial update. First improvement comes from multiplying the radial part 
of the Higgs field at all locations simultaneously by the same factor: -R(x) — > e^R{x.), 
where ^ is randomly chosen from a constant distribution around zero: ^ G [— e, e]. 
Under this update, the action ( |2.4| ) changes as AS'(^) = ae^^ + be^^ — a — b, where a 
and b are the sums of the terms proportional to R^ and R'^ in the action. Accounting 
for the measure factors, the update is accepted with the Metropohs probabihty p{^) = 
min(l, exp[4V^ — AS'(^)]), where V is the volume of the system. 
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The multiphcation factor has to be very close to unity in order for p(^) to be non- 
neghgible. In practice, we choose ^ from a window of width ~ 0.05 - 0.001 around zero, 
depending on the volume. The width of the window is chosen so that the acceptance 
is approximately 60-70%. 

In fig. we show the autocorrelation functions for an m*jj = 60, /^g = 8 system on a 
14^ lattice with both (a) local Metropolis and (b) local + global update. Even though 
the change in R{x) in each global multiplication is very small, the gain over only local 
Metropolis/heat bath is about a factor of 5. The additional cost of the global update 
in terms of cpu time is negligible, since it involves only one accept/reject step for the 
whole volume. 

Higgs field overrelaxation. Let us parametrize the Higgs field as $ = a^r^, /i = 
. . . 3, G R, where Tj are Pauli matrices and tq is the 2x2 unit matrix. From 
eq. (|2.4|) we see that the local potential of the Higgs field at location x is 

l^[$(x)] = -«^(x)F,,(x) + i?(x)2 + /3ij(i?(x)2 - 1)2 (3.1) 

where = a^a^ and is the sum of the hopping terms 

F^(x)=/?H^Trr^[ 5: $t(x - z)f/,(x - z) + $t(x + z)tA/(x)]. (3.2) 

^ j=l,2,3 

This form seems to suggest separate update steps for the radial and SU(2)-components 
of the Higgs field. However, this is not the optimal method: even though the overrelax- 
ation for the SU(2) direction can be readily performed, updating the radial component 
becomes quite complicated. In this case the overrelaxation step R ^ R' would require 
finding R' so that 

[dG{R)/dR]-^ exp[-G(i?)] , where G{R) = V{R) - logR^ , (3.3) 

remains invariant (here we use a notation where the x-dependence of the variables 
is suppressed) . This can be approximated by finding a solution R' to the equation 
G{R') = G{R), and performing a Metropolis accept/reject step using the probabil- 
ity weight [dG{R)/dR]/[dG{R')/dR']. This update has been used in 4d Higgs model 
simulations ||25|, However, the acceptance rate is only ~ 80%, and the algorithm 
behaves dynamically rather like a heat bath or Metropolis update. 

A more efficient method is to update the Higgs variables in the plane defined by 
4-dimensional vectors and F^, using the Cartesian components of parallel and 
perpendicular to F^: 

X = aJ,, Y^ = a,-Xf^, (3.4) 
where = F^/F and F = ^Jf\^. In terms of X and eq. ( ^31) becomes 

Vm = -XF + (1 + 2(3r{Y^ - 1))X^ + Pr{X^ + Y^) . (3.5) 
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Now we can update the X and Y components of $ separately: overrelaxation in Y is 
simply the reflection Y^ ^ Y^ = —Y^, or a'^ = —a^ + IX (this is exactly equivalent 
to the conventional SU(2) reflection overrelaxation procedure). In order to perform an 
exact overrelaxation to the X-component, we need to flnd X' so that 

[dV{X')/dX']~^ exp[-y(X')] = [dV{X)/dXY^ exp[-V(X)] . (3.6) 

To solve this equation we would have to resort to iterative numerical methods, which 
can be costly in terms of cpu-time. Instead, we used the following approximation to 
the overrelaxation (|3.6|): we flnd the solution to equation V{X') = V{X) and accept 
X' with the probability 

/ N dV(X)/dX , , 

V{X ) = min(po, 1) , Po = ^^F^^Tyy^ • (3-7) 

Since V[X) is a fourth order polynomial, solving the equation V^(X') = ViX^ boils 
down to flnding zeros to a third order polynomial (we already know one zero X' = X, 
which can be factored out). In all realistic cases the parameters of V{X) are such 
that there always is only one other real root, and it is straightforward to write a 
closed expression for X'. This update is an almost perfect overrelaxation: in our 
simulations the acceptance rate varies between 99.7% - 99.98%, depending on the fie 
used. The acceptance is high enough so that the "diffusive" update dynamics inherent 
in the Metropolis acceptance step does not play any role, and the evolution of the fleld 
conflgurations is almost deterministic. A different but related overrelaxation schema 



to the one described here has also been used in 4d simulations [221. 



Wavefront update. Let us next consider the order in which the lattice variables 
are traversed. In the conventional even-odd update, the lattice is decomposed into 
even and odd sectors - in our 3d case, to two sets with x + ?/ + 2; even and odd. The 
variables are flrst updated on all the even points before updating the odd points. Using 
this schema it takes L/2 whole lattice updates before any kind of signal from a given 
location on the lattice can propagate through the whole volume. 

In the wavefront traversal we pick an arbitrary 2-dimensional plane from the system, 
deflned by one of the conditions ±x ± y ± z = const, (modulo periodic boundary 
conditions). We select one of the directions perpendicular to the plane to be the positive 
direction. The update proceeds in two stages: flrst, we update the Higgs variables $ on 
this plane, and second, the gauge matrices f/ on the links emanating from this plane to 
the positive direction. Both updates are performed with the overrelaxation algorithms. 
After the gauge fleld update we move up to the next plane to the positive direction, and 
start the updates again. With this method, a planar wave of overrelaxation updates 
propagates through the volume in a single update sweep. In our implementation we 
keep the same orientation of the plane until the volume has been swept through 4-9 
times. 
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Figure 1: The autocorrelation function of the observable L = V\x.)Ui{x)V {'x + i) 
calculated from an m*fj = 60GeV, V = 14^ lattice, (a) Heat bath/Metropolis, (b) Heat 
bath/Metropolis with global i?-update, (c) 4x (overrelaxation with even-odd traversal) 
+ Ix (heat bath + global update), and (d) 4x (wavefront overrelaxation) + Ix (heat 
bath + global update). In all the cases, one update means going once through all the 
lattice points. 

The particular diagonal orientation of the plane is chosen in order to simplify the 
update: the spins and gauge links on the plane can be updated independently of each 
other, and all the gauge links have an equal footing with respect to the plane — there 
are no link variables within the plane, as would be the case if the plane was oriented 
along the principal axes. 

In fig. |T] we compare the autocorrelation functions of the standard even-odd overre- 
laxation (c) and the wavefront overrelaxation (d). In both cases we perform four overre- 
laxation sweeps through the volume, followed by one heat bath/Metropolis update and 
one global radial update. For both cases the individual Higgs variable overrelaxation 
step is the XY overrelaxation described in paragraph II above. The overrelaxation 
methods perform much better than the pure heat bath algorithms (a and b). The 
wavefront overrelaxation has much better initial decorrelation, as seen from the very 
rapid decrease in the autocorrelation function, even though it seems to have roughly 
the same exponential autocorrelation time as the conventional even-odd method. How- 
ever, this rapid initial decrease means that the integrated autocorrelation time r^^t is 
small, giving correspondingly small statistical errors for the observables. In our tests 
the wavefront overrelaxation had typically 1.5-3 times smaller Tint than the even-odd 
overrelaxation. 



m* = 60GeV ^^ = 8 V=14 
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Multicanonical update. Multicanonical update |^ is essential for the interface 



tension calculations with the histogram method (see Sec. p.3| ). At the transition tem- 
perature the probability distribution of some order parameter, say i?^, has two distinct 
peaks corresponding to the two pure phases. The probability that the system resides 
in a mixed state, consisting of domains of the two phases separated by phase interfaces, 
is suppressed by the interface tension times the area of the interface (see, for example, 
fig. I). 

To enhance the probability of the mixed states the action eq. ( p. 41) is modified with 
the multicanonical weight function W: 

Smc = S + W{Rl) , = E ^ W (3-8) 

X 

The weight function iy(_R|) is chosen so that the resulting distribution p{Rs) is ap- 
proximately constant in the interval < i?| < R^2^ where and R^2 denote 
the pure phase peak locations. This is the main disadvantage of the multicanonical 
method: a priori, the weight function is not known; an exact knowledge of the weight 
function is equivalent to knowing the probability distribution of the order parameter, 
which is one of the quantities we attempt to calculate with the Monte Carlo simulation. 

The canonical expectation value of an operator O can be calculated by reweighting 
the individual multicanonical measurements Ok with the weight function: 

where the sums go over all measurements of O and 

The choice of i?| for the argument of the weight function is by no means unique; 
equally well one could use, for example, the hopping term J2x,i 2"^^ $^(x)f/j(x)$(x + i). 
The advantage of i?| is that in this case the weight function does not modify the update 
of the gauge fields U and the SU(2)-direction of the Higgs field V. However, the XY 
overrelaxation described above has to be modified. 

We parametrize W with a continuous piecewise linear function: 

WiRl) =w^ + (w,+i - Wi) ^^s-n ^ ri<Rl<n+i. (3.10) 

ri+i - Ti 

An initial guess for the parameters Wi can be obtained by finite size scaling the prob- 
ability distributions obtained from simulations using smaller lattice sizes. If deemed 
necessary, the parameters are further adjusted after preliminary runs. 

Since we perform most of the simulations with vector supercomputers, it is important 
to vectorize the multicanonical update. This is achieved with the following steps: 

(i) If Tj < i?| < Tj+i initially, the weight function is fixed to the linear form 
W'{Rl) = Wi + {wi+i - Wi){Rl - ri)/(ri+i - r^) for all 



9 



(ii) A number of Higgs variables are updated with the action S' = S + W . Since W 
is hnear in E?, it is straightforward to do this with a fully vectorized algorithm. 
We update $ in ~100-200 points in one vector. 

(iii) The whole vector of updates is accepted with the probability = min(l, 
eMW'{Rl)-W{Rl)]). 

Obviously, the acceptance in the step (iii) decreases when the number of variables 
updated in a single vector increases. Here the vector length ~ 100 is short enough so 
that i?| in practice often remains between the limits and rj+i after the update, and 
the acceptance is exceptionally good: in our runs the rejection rate was only ~ 10~^! 
At the same time the vector length is long enough so that increasing it does not give 
any significant gain in computational speed. Thus, the performance hit caused by the 
implementation of the multicanonical algorithm is negligible. 

The simulations were performed with Cray C-90 and X-MP supercomputers and, 
for the smaller volumes, with IBM RS6000 and HP 9000/735 workstations. The total 
amount of computing power used was about 5 ■ 10^^ flop = 160 Mflops year. 



4 Properties of the broken phase: the 3-loop effec- 
tive potential 

We shall first study the 3d SU(2)+Higgs system when it is in the broken phase. This 
will permit us to show the accuracy of the method and, in particular, to determine 
the size of the so far uncomputed 3-loop term in the effective potential (see (1.73), 
(11.151-152)), and verify that it is linear in 0. 

Although the object of study is the 3d theory, we use the simple parametrisation 
in eqs. ( p.8| - p.l(i| ) to permit one to use the 4d quantities m|^, T* in fixing the 3d pa- 



rameters x,y. Whenever m*fj,T* are quoted, the values of gi,x,y are to be computed 



from eqs. (|2.8| - p.l0|) . The same values of m*jj,T* can correspond to different physical 
parameters in different theories. 

We study the effective potential in the broken phase with the help of the scalar 
condensate H- Eq. (11.139) relates the value of the scalar condensate in the MS 

scheme, where 

(0t0(;,,)) - (0V(/i2)) = ^ log ^, (4.1) 

lOTT^ /i2 

to the corresponding lattice quantity (-R^). The relation, in which the correction term 
vanishes in the continuum limit, is 



« . 1 ((« V 4) - (log + C + - ^) + o(-l 

^1 8 V tcPhJ (47r)2V 2^ 4 J 
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110 


145 


155 


165 


167 


{(t)^{T*))/T* 


1.9553(13) 


0.7061(7) 


0.4718(6) 


0.2506(10) 


0.2052(36) 


jj, = 1.2mT 
jj, = A.TrriT 


1.9642 
1.9650 


0.7184 
0.7207 


0.4861 
0.4898 


0.2707 
0.2793 


0.2235 
0.2359 


CW/i = 1.2mT 
CW/i = 4.7mT 


1.9649 
1.9660 


0.7206 
0.7223 


0.4902 
0.4922 


0.2818 
0.2847 


0.2397 
0.2433 



Table 1: Measured and computed values of (0V(T*))/T* for = 80 GeV and 
T* = 110, 145, 155, 165, 167 GeV. The first row gives the lattice value in the limit 
V oo, a — s> 0, obtained using eq. ( [4.21 ). The following two rows with /i- values refer to 
the result in h expansion including the 2-loop term and the known part of the 3-loop 
term (Appendix A, eq. ( |A.22| ) with P{h) = 0), and the last two rows to a numerical 
CW calculation using the RG-improved 2-loop effective potential. 



Here numerically ( + ~ S = 0.67. Thus, the extrapolation of lattice measurements 
of the quantity (i?^) to the limit jSc ^ oo allows one to determine an "exact" value of 
the scalar condensate which then may be confronted with the perturbative expression. 

We use the strategy appropriate for simulations at T < Tc as explained in Sec. 0. 
We take Tif^ = 80 GeV and select several temperatures for which the system is in 
the broken phase: T* = 110, 145, 155, 165, 166, 167 GeV. For each temperature we 
choose different values of Pg (usually, 6 < /5g < 32). These numbers completely define 
the parameters of simulation through the constant physics curve ( |2.5| - |2.7|) . For every 
simulation, the volume of the system was large enough to make finite size effects smaller 
than the statistical uncertainty in (-R^). 

Examples of measurements at T* = 110, 145 GeV and 165 GeV are shown in figs. ^ 
^, and the values of {(f)^(f){T*)) /T* in the continuum limit are collected in table 1. The 
range of the temperatures used is quite wide, the expectation value of the Higgs field 
varies (in 4d units) from 0b/T* ^ 2 to 4>b/T* ^ 0.6. 

The perturbative 2-loop computations were done with two methods, described in 
Sec. 5 of II. The first method is based on a straightforward h expansion of the conden- 
sate. The result is given in Appendix A (eq. (|A.22|) ) and shown in fig. Fig. ^ contains 
the tree, 1-loop and 2-loop terms (eqs. (11.76-78)) but also that part of the 3-loop term 
which is known because it is related to lower order potentials. In the second method 
(which we call the CW - Coleman- Weinberg - type of computation), one numerically 
finds the location of the broken minimum of the 2-loop effective potential, determines 
the ground state energy, and then computes the condensate with the help of (11.42) as 
dV{(j)b)/dm1. To avoid infrared divergences, the straightforward CW- method has to 
be somewhat improved, see [Appendix B . 

The expansion parameter of the h expansion is ~ S's/v— "^i and thus it converges 



the better the deeper one is in the broken phase. This also means that the signal for 
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Figure 2: Data for {(f)^(f){T*)) /T* as a function of 1/(3g for = 80, T* = 110 
GeV computed from measured values of (i?^) using eq. ( |4.2|) with = T*, = 
0.44015T* (parametrisation in eq. ( |2.8| )). The perturbative values corresponding to 
(3 = 0,20,40,50, calculated with the CW-method at the scale /i = 2.37mT and then 
run to n = T* with eq. (|4.1|), are shown on the vertical axis. 




Figure 3: As fig. |, but for T* = 145 GeV and T* = 165 GeV. 
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Figure 4: The values of {(j)^(j){T*)) /T* in the h expansion as a function of temperature 
for m*fj = 80 GeV from eq. (|A.22|) . The 3-loop curve contains the known part thereof, 
with l3{h) = 0. It is seen that the ^-calculation becomes increasingly unreliable as one 
approaches the critical temperature. 



the 3-loop term vanishes if one goes too deep into the broken phase. Since one cannot 
be too close to m| = (or T* = 173.5 GeV in table 1) either, there is an optimal region 
in between. The numerical CW method works also at Tc and is thus more accurate. 
We include the h results mainly since the formulas are very explicit. 

When all orders of perturbation theory are summed, {(j)^ (f>{fi)) depends on fi ac- 
cording to eq. (|4.1| ). At a finite order in perturbation theory, however, there is extra 
/i-dependence which can be used as an indication of the accuracy of the calculation. 
In table |l] the parameter /i used in the calculation of "^^is varied within the 

limits 0.5-2.37 < [x/rriT < 2-2.37, the central value 2.37 (table 2 of I) being determined 
by the requirement of the "best" convergence of perturbation theory for the effective 
potential in the vicinity of the minimum. Then (0^0(/i)) was run to fi = T* according 
to the exact running in eq. ( [4.1|) . 

Inspection of table |T] shows that the CW values of the condensate are practically 
independent of /i so that the convergence of CW type perturbation theory is better 
than that of the ordinary one. The difference between these two methods is due to 
higher order corrections; the CW method sums a subset of them and can be used at 
the phase transition, where the first method fails^. One can see that the CW values of 
the condensate are larger than the corresponding lattice values, and that the difference 



^The CW method provides an automatic summation of one-particle reducible diagrams, whereas in 
ordinary perturbation theory they have to explicitly calculated to the desired order, see Appendix A. 
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110 


145 


155 


165 


167 


Pew 


53.0(7.3)(2.7) 
58.0(7.4)(2.3) 


41.8(2.1)(5.1) 
48.1(2.3)(2.3) 


38.3(1. 3)(6. 7) 
46.6(1. 8)(3.6) 


35.1(1.1)(12.0) 
49.1(1.3)(2.6) 


28.8(2.8)(14.4) 
46.2(3.6)(2.5) 



Table 2: The values of the 3-loop coefficient /? at different temperatures using the 
h expansion and the numerical CW computation. The ffist number in brackets is 
statistical uncertainty and the second is an estimate of the systematic error associated 
with the change of the parameter /i within the interval in table ^ 



between the two is statistically significant. This proves the necessity of higher order 
effects and singles out their sign: the "exact" vev of the Higgs field at some fixed 
temperature is smaller than that given by 2-loop perturbation theory. 

The study of the difference between the 2-loop and lattice values of the scalar con- 
densate allows one to determine the magnitude and the structure of the 3-loop term. 
In CW perturbation theory for the scalar condensate we add to the 2-loop effective 
potential a linear term, expected on the 3-loop level (see (1.73)):^ 

m = j^,9tmT{<P). (4.3) 

Here /? is a constant to be determined by fitting the prediction to the lattice number 
at each temperature (for constant ^«|^)[]. In the h expansion method the procedure is 
extremely simple: one takes (0V(/^)) from eq. ( |A.22| ) of Appendix A, runs it to the 
scale T* using eq.( [4.1|) and fits the constant P to get agreement with data. In other 
words, (3 is linearly proportional to the difference of the perturbative and the lattice 
value in table |l|. The outcome of this procedure is in table |. 

Within the T* range discussed the vev of the scalar field varies from 0.6T* to 2T*. 
Since the value of (3 does not depend on T* within error bars, one may conclude on 
the basis of the lattice data that for around 80 GeV there exists in the effective 
potential a 3-loop linear term with positive sign and /5 ~ 50. 

We estimated /? at different Higgs masses also by another method. We computed 
on the lattice the continuum limit of the critical temperature and the value of the 
condensate at T* (the methods are described in detail in the next section). We found 
for m*jj = 60 GeV that T* = 138.38(5) GeV and {(p^(p{T*)) /T* = 0.227(6); and for 
m*H = 70 GeV that T* = 154.52(10) GeV and {(ji^^T*)) /T* = 0.162(12). These 
numbers may be reproduced with the 3-loop effective potential with (3 ~ 49(2) and 
with (3 ~ 46(4), respectively. The small Higgs mass = 35 GeV is not informative 

^ The structure of the 3-loop effective potential is discussed in more detail in Appendix A. 

^In Q a rough estimate of the parameter /3 was given, /3 = —15(20) with only statistical errors 
quoted. After the analytical computation of the constant physics curve the systematic errors can 
be removed, while the higher statistics allows to reduce the statistical errors considerably. 
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since higher order corrections at the critical temperature are numerically small and the 
extraction of (3 with any reasonable accuracy is not possible with the data we have. 

Results at m*^ = 60, 70 GeV indicate that within errorbars, P does not depend 
on m*jj. This is what one expects: in the abelian U(l)+Higgs model there is no 
linear term so that its coefficient should be proportional to the non-abelian gauge 
coupling. 



5 The phase transition 

Since we are mostly interested in the properties of the phase transition, most of our 
simulations are performed at and immediately around the transition temperature. We 
study three different Higgs mass parameters rriff = 35, 60 and 70 GeV. We have mostly 
concentrated on 60 GeV Higgs, since it is close to the physically allowed mass range, 
but the transition is not yet too weakly ffist order to be studied with moderately sized 
lattices. 
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Table 3: Lattice sizes used for the simulations at the transition temperature for each 
(mJ^,/3G)-pair. In most of the cases, several /^//-values were used around the transition 
point. Multicanonical simulations are marked by subscript (m)- 

For each of the three values of we use the gauge couplings Pg = 8, 12 and 
20, and for mj^ = 60 GeV also Pg = 5 (remember that Pg is directly related to the 
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lattice spacing through Pg = 4/((7|a)). In table || we list the lattice sizes for each 
{m*fj,f3G) -pair. Each lattice has several runs with different values for Ph in order to 
accurately locate the transition; typically smaller lattices have 5-20 and larger ones 
1-3 values of /3j^ . Separate runs are then joined together with the Ferrenberg-Swendsen 
multihistogram method . When the interface tension causes noticeable supercritical 
slowing down, the multicanonical algorithm is used. All in all, the total number of 
separate 'runs' — different combinations of lattice sizes and coupling constants — 
described in this paper is 289; this includes also runs away from the transition region. 

The large number of lattice volumes with several lattice spacings makes it possible 
to accurately extract continuum values of the physical observables: 

(1) For fixed lattice spacing a (fixed we extrapolate the lattice measurements 
to the thermodynamical limit V —>■ oo. 

(2) Each of the V = oo values are in turn extrapolated to the continuum limit a ^ 
(Pg - oo). 

Note that in the 3d continuum theory the lowest-dimensional gauge-invariant opera- 
tor 0^0 has the dimensionality GeV. This means that the scaling violations in physical 
quantities start from the first power of the lattice spacing a (in the 4d theory, scaling 
violations are proportional to a^). 

The transition becomes weaker — the latent heat and the interface tension become 
smaller — when increases. Measured in dimensionless lattice units, the transition 
also becomes weaker when the lattice spacing a decreases for fixed m}j. This can be 
observed from the probability distributions of the average Higgs field squared: = 
yJ2xR(?^)'^- In fig- H we show the distributions for some of the largest volumes for 
m*jj = 60 GeV, Pg = 5, 8, 12 and 20. When Pg increases, the separation between the 
peaks becomes smaller and the minimum between the peaks becomes shallower. As we 
will explain in Sees. |5]^ and ^73| , these features are directly related to the latent heat 
and the interface tension, respectively. 

In fig. P the corresponding Pg = 8 histograms are shown for = 35 and 70 GeV. 
The dramatic effect of the Higgs mass to the strength of the phase transition is clearly 
evident. Note that all the distributions in figs. |^ and |^ correspond to different values 
of Ph', each of the histograms has been reweighted to the "equal weight" value (see 
Sec. P). 



5.1 The critical temperature 

The critical temperature can be determined extremely accurately from the Monte Carlo 
data. In the SU(2)+Higgs model there are no known local order parameters, which 
would acquire a non-zero value only in one of the two phases of the model. Instead, 
we use order parameter like quantities which display a discontinuity at the transition 
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Figure 5: The probability distribution of the average Higgs length squared for 
m*H = 60 and (3g = 5, 8, 12 and 20. 



point (when V ^ oo). The quantities we use are and the hopping term 

L 



^$:^Tr0(x)t/.(x)y(x + 



(5.1) 



where V^(x) is the SU(2) direction of the Higgs variable $(x) = R{x)V{x.). The 
behaviour of (L) as a function of Ph is shown in fig. Ul for = 60 GeV, /3g = 8 case 
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Figure 6: distributions for [3g = 8, m*^ = 35 and 70 GeV. 



for lattice sizes up to 32'^. The development of the discontinuity is clearly visible. The 
continuous lines are a result of the multihistogram method calculation. 

For each individual lattice volume and (m^, /9G')-pair, we locate the pseudocritical 
coupling (3h,c with several different methods (see, for example, ||2^ and references 
therein): 

(1) maximum of the L-susceptibility C{L) = {{L — {L))"^) 

(2) maximum of C{R^) = {{R^ - {R^)f) 

(3) minimum of the 4th order Binder cumulant of L: B{L) 

(4) "equal weight" Ph -value of the distribution p{R^) 

(5) "equal height" Ph -value of the distribution p{L) 



1 - {L')/{3{Lr' 



The locations of the extrema of the observables C{L), C{R^) and B{L) are computed 
by reweighting the original measurements; the error analysis is performed with the 
jackknife method, using independent reweighting for each of the jackknife blocks. As 
an example, we show C{R^) and B{L) for m*^ = 60 GeV, /?g = 8 lattices in fig. ^ as 
functions of (3h- For clarity, the errors are omitted from the figures. Again, we would 
like to point out the unambiguous first order scaling displayed by the data in these 
figures. 

The /^H-values for the "equal weight" and the "equal height" distributions are also 
found by reweighting the histograms independently for each jackknife block. The his- 
tograms in figs. 1^ and ^ are all equal weight R^ histograms; that is, the areas in the 
symmetric phase and broken phase peaks are equal. The calculation of the peak area 
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Figure 7: The hopping term (L) as a function of (in for different lattice sizes, for 
= 60 GeV, /5g = 8. 

requires an arbitrary selection of the value of which is used to separate the peaks; 
we used a fixed value for all the lattices in each (m^, /3G)-set, defined by the minimum 
of the distribution of the largest volume. 



The infinite volume limit: The values of ^h,c determined with the methods (l)-(5) 
above differ for each individual lattice, but the V ^ oo extrapolations are very well 
compatible within the statistical errors. It should be noted that the different methods 
for determining fiji^f. do not give statistically independent results, and it is not justified 
to combine the values given by different methods together. However, they serve the 
purpose of checking the consistency of the infinite volume and continuum limits. 

In figs. Pl-pm we show the infinite volume limits for mj^ = 60 GeV and /3g = 8, 12 
and 20. As can be observed, different methods converge extremely well (the intercepts 
of the dashed lines at XjV = are nearly equal). The same holds true for Higgs 
masses mj^ = 35 and 70 GeV, which are not shown here. In table ^ we show the 
V CO extrapolations of Ph,c, using the data obtained with the equal weight of p{R^) 
-method. The corresponding values of the critical temperature T^*, calculated with 
eqs. ' 



( |2.10|) , are also shown. 



The continuum limit: In fig. 11 the V = oo values of T* are extrapolated to the 



continuum limit a ^ for m*^ = 60 GeV. We expect the leading deviation from the 
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Figure 8: C(i?^), the susceptibility of the order parameter i?^ (top), and the Binder 
cumulant B{L) of the order parameter L (bottom), plotted as a function of Ph for 
m*fj — 60 GeV and /3g = 8 runs. 



m^/GeV 


Pg 




T*/GeV 
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20 


0.3391279(26) 


154.03(16) 



Table 4: The infinite volume critical couplings (3h,c, and the critical temperatures T*. 
The values of (3h,c are calculated from the "equal weight of p(i?^)" data. 

continuum limit value to be of order C(a); in this case, there are 4 values for Pc, and 
the accuracy of the data is high enough that quadratic fits are needed in order to have 
good x^/d.o.f. values for the fits. The quality of the fits is very good for the T* values 
calculated with any of the five criteria, and the final extrapolations are statistically 
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Figure 9: The V = oo Umit of the pseudocritical coupUngs (3h,c for m*^ = 60 GeV, 
(3g = 8, calculated with five different methods. All the methods give compatible 
= oo limits; the point near 1/V = is the result from linear extrapolation of the 
"equal weight of p^R"^)" -values. Only the largest volumes are shown. 

compatible. For concreteness, we use the equal weight of p{R^) -results for our final 
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Figure 10: The same as fig. ^ for /^^ = 12 and /^g = 20. 
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60 GeV. A quadratic fit (dotted line) in 1/Pg to all four values gives good x^/d.o.f., 
whereas a linear fit (dashed line) using only /3g > 8 does not have acceptable x^- The 
extrapolations shown are for "equal weight of p{R^y^ data. 



numbers. 

In fig. |l^ we show the corresponding extrapolations for nf^ = 35 and 70 GeV. In 
these cases we use linear fits. The final results are summarized in table ^ together 
with the perturbative values of T*. 





T;/GeV 


T;P<=''VGeV 


35 


92.64(7) 


93.3 


60 


138.38(5) 


140.3 


70 


154.52(10) 


157.2 



Table 5: The continuum limit extrapolations of the critical temperatures. The 
60 GeV point has been calculated with a quadratic fit, others with linear fits. 



There is a systematic difference between the perturbative and the lattice results; 
the perturbative T* is considerably larger. The values are closest to each other when 
m*fj = 35 GeV, but for m*j^ = 60 and 70 GeV the difference is more than 20 standard 
deviations. These results agree qualitatively with the results from the simulations 
with the Ao field |l|, §, but the errors here are almost an order of magnitude smaller. 
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Figure 12: The same as fig. 11 for m*^ = 35 GeV and 70 GeV, using linear fits. 



They also agree qualitatively with 4d simulations ^ and recent 3d simulations by 
Ilgenfritz et al. [|15|. Note, however, than in 4d case the errors are considerably larger. 



and in [|1^] the extrapolation to continuum limit is not taken. 



5.2 The latent heat and v{T*] 
The latent heat L — the energy released in the transition — can be calculated from 



L dAp T d ^ ^ ^ T d ^ ^ 

— = — - = A bg Z = AP . 

T dT VdT ^ VdT 



(5.2) 



where the derivatives are evaluated at the critical temperature, Ap is the difference 
of the pressures of the symmetric and broken phases, and AP is the difference of the 
probabilities of the phases in volume V. In eq. ( |5.2| ), T is the physical (4d) temperature; 
for simplicity, in the following we substitute T T*. In Sec. we discuss how the 
correct physical result can be obtained. The quantity AP is directly proportional to 
the difference of the areas of the two peaks in the order parameter distributions near 
T*, and d{AP)/dT* is readily calculable by reweighting. 

An alternative method is to evaluate d{Ap)/dT* directly from the action in eq. ( |1.1|) : 



L 

c 



m 



T, 



*3 



A( 



1 



^j;§9ipHf3GA{R') 



(5.3) 



Both methods give compatible results in the V 
have been calculated with eq. 



oo limit. The results shown here 
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Figure 13: Extrapolation of the latent heat to ^ oo limit (left figure) and to the 
continuum limit a — (right figure) for = 60 GeV lattices. 



Fig. displays the limits V ^ oo and a for the latent heat for = 60 GeV 
systems. The infinite volume limit is taken by extrapolating linearly with respect to the 
inverse area 1/(AT*^) = (L^^aT*)"^ of the system, where L^a is the linear length of the 
lattice (one of the short dimensions for the cylindrical volumes). The extrapolation with 
respect to the inverse volume would fail to accommodate the cubical and cylindrical 
lattices simultaneously with the same scaling ansatz. The inverse area -type behaviour 



of the latent heat is known to occur for the Potts models in 2 dimensions pO| , pT] (the 
area in this case being the linear dimension of the lattice). 

The Higgs field expectation value in the broken phase v{T*), defined here by 

(5.4) 



can be calculated from (i?^) using eq. ( [4 .21) . Because of the close relation between the 



equations ( |5.2| ) and ( [4.2| ), the limits are very similar in both cases. In fig. |T^ we show 
the approach of v'^{T*) to the continuum limit. 

Higgs masses 35 and 70 GeV are analyzed in a similar way; the results for v{T*) and 
the latent heat are shown in table ^ 

Note that the values for both v and L are quite close to the perturbative values, 
evaluated at the perturbative critical temperature T*^. If we use the lattice critical tem- 
perature f P becomes larger, as can be seen from the last column in table This shows 
the presence of higher-loop perturbative corrections in the broken phase, discussed in 
detail in the previous section. 
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Figure 14: The continuum limit of the square of the Higgs field vev f^(T*). 



m^/GeV 


l/{t:y 


Lp/(T;p)^ 


v{t:)/t: 




v^{t:)/t: 


35 


0.256(8) 


0.22 


1.86(3) 


1.75 


1.87 


60 


0.0406(7) 


0.040 


0.674(8) 


0.68 


0.82 


70 


0.0273(16) 


0.027 


0.57(2) 


0.55 


0.70 



Table 6: The latent heat L and the Higgs field expectation value v{T*) in the broken 
phase. Here T*^ is the perturbative critical temperature (see table ||). The comparison 
between lattice results and perturbation theory is discussed in Sec. 0. 



5.3 The interface tension 

The interface tension is one of the primary quantities which characterize the strength 
of the phase transition. We measure it with the histogram method [^: at the pseud- 
ocritical temperature, a system in a finite volume predominantly resides in either the 
broken or the symmetric phase, but it can also exist in a mixed state consisting of do- 
mains of the two states. The probability of the mixed state is suppressed by the extra 
free energy associated with the interfaces between the phases. This causes the typical 
2-peak structure of the probability distribution of the order parameter at the critical 
temperature (see figs. | and H): the midpoint between the peaks corresponds to a state 
which consists of equal volumes of the symmetric and broken phases. Because of the 
associated extra free energy, the area of the interfaces tends to minimize. Assuming 
a lattice with periodic boundary conditions and geometry x L^, where < Lz, 
the minimum area is 2 x A = 2(L^a)^ — the number 2 appears because there are two 
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Figure 15: The interface tension for = 60 GeV extrapolated to ^ cxo (left) and 
a — i> (right). 



separate interfaces. The interface tension cr can be extracted from the limit 



a 
f 



lim 



1 



v^oo 2A 



log 



Prr 



(5.5) 



where Pmax and Pmin are the probability distribution maximum and the minimum 
between the peaks. 

In practice, the infinite volume value of a is reached in such large volumes that careful 
finite size analysis of eq. ( ^.5|) is necessary. Numerous studies exist in the literature 
33, H, p5l; here we follow [0: 



a- 



T 



log 



1 



max 
^ mm 



log -^2 - ^ log Lx + ]^G + const. 



(5.6) 



The function G interpolates between lattice geometries; the limiting values are G = 
log 3 for cubical volumes (L^ = L^) and G = for long cylinders (L^ ^ L^). 

The finite size scaling ansatz ( |5.6| ) assumes that the two interfaces are far enough 
apart from each other that their mutual interaction is negligible. In practice, this is very 
difficult to achieve in cubical volumes and usually requires the use of long cylindrical 
lattices. The order parameter histograms develop a flat minimum when this condition 
is fulfilled: the flatness signals a constant free energy when the volume fractions of the 
two phases slightly change, and the interfaces move with respect to each other. 

In order to find the fiat part of the histograms we use large cylindrical volumes for 
all jSc values for = 60 GeV. The volumes included for the analysis are: jSc = 5: 
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122 X 72 and IQ^ x 80; /3g = 8: 20^ x 140, 24^ x 120, and 30^ x 120; /3g = 12: 26^ x 156, 
30^ X 150, and 36^ x 144; (3g = 20: 40^ x 200 and 50^ x 200. In fig. | we show the 
equal weight histograms of ^(-R^) for these lattices. 

In the interface tension analysis we use equal height histograms of p{L) . Equal 
height histograms were chosen instead of equal weight because the determination of 
Pmax in equation ^ then becomes unambiguous; and p{L) instead of p{R^) since the 
former histograms have more symmetric shapes than the latter, and the equal height 
and equal weight I3h,c values are quite close to each other. Nevertheless, we check 
the measurements with equal weight p{R^) histograms (using Pmax which is a linear 
interpolation of the two peak heights to the B? value of the minimum); the results are 
very well compatible within statistical errors. 

A comment about extracting the extrema from the histograms is in order: we find 
the maximum values by fitting a parabola close around the peaks of the histograms, 
and the value of the fiat minimum by fitting a constant. This method gives much 
smaller errors and more reliable results than simply using the absolute extrema values, 
which are very prone to statistical noise. The same method was used also to locate the 
equal height /Sjy^c-values. 

In the left part of fig. 



15 we show the interface tension measurements from each 



lattice. The values shown here include the finite size scaling correction l/L^(|logL2 — 
llogLj.) from eq. (p.6|), so that only a factor const./ Li remains. For each /3g, the 
behaviour of the data is linear in 1/L^. On the right part of fig. 15 we extrapolate a 
to the continuum limit. The results are shown in table |^. 





a/{T:f 


^p/(T;p)3 


35 


0.0917(25) 


0.066 


60 


0.0023(5) 


0.0078 


70 




0.0049 



Table 7: The interface tension a. Only the m*fj = 60 GeV result is an extrapolation to 

35 GeV value is only from (3g = 8 simulations. 



the continuum limit; the m*jj 



The m}j = 60 GeV continuum limit result in table |^ is obtained by a linear extrap- 
olation in I/Pg- However, on closer inspection the data from Pg = 12 and 20 lattices 
on the left part of fig. |TB| seem to indicate that the interface tension is already scaling 
when Pg > 12, even in finite volumes. If we use the finite volume (3g = 8 and 12 data 
and extrapolate to infinite volume, we obtain the result a/{T*)^ = 0.0042(3), which is 
not compatible with the result in table 0. Nevertheless, since the data for the latent 
heat or v{T*) do not display similar scaling behaviour, we use the linear extrapolation 
in I/Pg in si'll cases. 

For m*jj = 35 GeV, only (3g = 8 lattices are cylindrical enough so that we can 
estimate the V ^ oo limit. However, the continuum value can not be extrapolated. 
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For m*jj = 70 GeV we did not observe good enough flat parts in the order parameter 
histograms, and we do not attempt to estimate the interface tension. 

6 The correlation lengths for m*^ = 60 GeV 

The measurement of the Higgs and W masses around the transition temperature is 
somewhat orthogonal to the measurements of the quantities directly associated with 
the transition itself: instead of attempting to enhance the tunnelling of the system 
from one phase to another, in this case it is imperative that the system stays in one 
homogenous phase throughout the measurement. This is due to the spurious signal 
caused by the tunnelling correlations: mtunnei oc exp(— cxA/T) |2^. Even an incom- 
plete tunnelling can make the extraction of the physical mass very complicated. Since 
the tunnellings are suppressed by the exponential factor exp(— crA/T), we perform a 
separate set of simulations around the critical temperature using large volumes and 
monitor the simulation time history of order parameters in order to ensure that the 
system stays in a single phase throughout the measurement. 

We perform the measurement only for = 60 GeV systems, using lattice sizes 
30^ X 60 and 40^ x 80 for (3g = 8 and 40^ x 80 and 50^ for Pg = 12. Let us define 
operators 

Cd(^)= E E^:^'3..Tr[r'^$^(x)f/,(x)...?7,(x+(ci-l)e,)<l>(x + cie,)] (6.1) 

i=l,2 X 

where a = 0, ... 3; r", a = 1, 2, 3, are the Pauh matrices and r° = 1; and ci = 0, . . . , 4 
is the length of the 'chain' of link matrices between $^ and $. Using eq. (|6.1|) we can 
define correlation functions sensitive to the Higgs and W channels: 

h,{l) = + (6.2) 

Ml) = T^EEcS(^K(^ + 0- (6.3) 

z a=l 

The masses mniT*) and mw[T*) are found from the exponential fall-off of hd and 
Wd- The results are independent of the parameter d, and it is chosen to minimize the 
statistical errors. We obtain best results with d = for hd (c° = $"1"$) and d = 4 
for Wd- 

The measured values of mH(T*) and mw(T*) are shown in fig. The scaling 
between Pg = 8 and 12 is very good, and both mn and mw display a discontinuity 
at the transition. At the transition temperature we are able to maintain the systems 
in either the broken or the symmetric phase throughout the measurement, so that 
immediately around we have two values for the masses. Both m/^(T*) and m-\^[T*) 
are higher in the symmetric phase than in the broken phase. 
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Figure 16: The Higgs masses (top) and the W masses (bottom) for m*jj- = 60 GeV, 
/3g = 8 and 12 systems as functions of the temperature. The perturbative results are 
shown for rriH with sohd (tree-level value), dashed (1-loop pole) and dot-dashed (some 
2-loop corrections, see Sec. |10.1| ) lines; for mw with solid (tree-level) and dashed (1- 
loop) lines. Note that the ratio mw/rriH is much larger than indicated by the tree- level 
curves, since radiative corrections make mw larger and tjih smaller. 
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7 The metastability ranges and scaling properties 
of the symmetric phase 

An important characteristic of the phase transition is the range of metastabihty of 
the broken and symmetric phases. The symmetric phase is metastable below the 
critical temperature down to T = T_, and the broken phase is metastable above the 
critical temperature up to T = T+. The aim of the present section is to estimate 
the metastability range, and to study the dynamics of the composite field by 
determining its effective action. We shall concentrate on the case m*fj = 60 GeV 
(x=0.06444) and T* near T* {\y - ?/,| < 0.06). 



7.1 The metastability range from the correlation lengths 



Consider the behaviour of the scalar mass near the phase transition, fig. |T6|. When 
the temperature decreases, mn in the symmetric phase rapidly decreases. This be- 
haviour suggests that it reaches zero at some point, and the symmetric correlation 
length diverges. This is the lower spinoidal decomposition point and corresponds to 
the temperature at which the symmetric phase ceases to be metastable, T_. Similar 
behaviour takes place in the broken phase, when the temperature increases towards T+. 
This observation allows us to estimate the metastability range, i.e., the upper and lower 
spinoidal decomposition temperatures. In general, one would expect the following de- 
pendence of the Higgs correlation lengths on the temperature: 

CH = Q\T*/T,~ir\ (7.1) 

where i = b,s labels the broken and symmetric phases, and Tb = T+,Ts = T^. The data 
we have does not allow one to determine the critical exponents and the temperatures 
T± simultaneously with good accuracy, so that we have chosen 7j = 1/2, following the 
guidance from mean field theory. The results for the metastability region are given in 
table H, using the = 60 GeV, /3g = 8 and 12 results from fig. |16[ 



Pg 


T_/GeV Csgi 


T+/GeV c^gl 


8 
12 


135.0(4) 0.35(2) 
135.5(3) 0.33(2) 


140.1(4) 0.49(2) 
139.6(4) 0.45(2) 



Table 8: The endpoints of the metastability temperature ranges for the symmetric (T_) 
and broken (T+) phases, determined from = 60 GeV Higgs correlations with the 
ansatz (fTll) . 
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Figure 17: Evolution of the potential t/(vr) for the effective field vr = related to 
the lattice variable B? through relation ( [4 .21) , in the vicinity of the phase transition for 
m*^ = 60 GeV and /3g = 8. Different Ph values are reached by reweighting a single 
(multicanonical) simulation. 



7.2 Metastability from reweighting 

A direct method to measure the metastability range is to use the order parameter 
distributions and reweighting. Let us define an effective variable tt = = ^R'^ + 
const., where the last equality follows from eq. ( [4.2p . Around the pure phase peaks, 
the probability distribution p{R^) is related to an effective potential U{n): 

piR") oc e-^^W , (7.2) 



where V is the volume of the system. Eq. (|7.2| ) has pre-exponential corrections; how- 
ever, to the accuracy we are working here the above formula is sufficient PBI. By 



reweighting the distribution p{R^) we obtain the temperature dependence of the po- 
tential. It should be noted that eq. ( [7.2|) is only valid in the immediate neighbourhood 
of the pure phases; it does not correctly describe the mixed state between the pure 
phases. 

The endpoints of the metastability branches can be found by locating the temper- 
ature at which the barrier against the tunnelling vanishes — the minimum of the 
potential in eq. ( |7.2| ) turns into an inflection point. We present the evolution of the 



31 



139.5 
139.0 h 
138.5 
138.0 F 

> 

O 137.5 

* 

137.0 - 
136.5 
136.0 h 



135.5 



m* = 60 GeV 



-o- 



0.00 



0.05 



0.10 



0.15 



0.20 



Figure 18: The metastability range for = 60 GeV Higgs systems, calculated with 
the reweighting analysis and extrapolated to a — 0. 



potential in fig. 1/ for the symmetric and the Higgs phases for an = 60 GeV, 
(3g = 8 system. These potentials are reweighted from the 30^ x 120 lattice histograms 
shown in fig. ^ In fig. |1^ we extrapolate the temperatures T_ and T+ to the continuum 



limit, with the result 



Tl = 137.37(12) GeV T; = 138.72(15) GeV. 



(7.3) 



The metastability range here is smaller than the range determined by fitting the cor- 
relation lengths (table We believe that the values in eq. ( |7r3| ) are more rehable, 
because the power law extrapolation of the correlation lengths is very sensitive to sta- 
tistical errors and finite volume effects. Note also that the range of Ph over which 
reweighting is carried out is very small. 



7.3 The effective theory description 

More information on the dynamics of the composite field vr is contained in its effective 
action. We shall write it separately for each phase in the form 



cS 



(fx 



1 1 

Z22 



{d.nf + U{n) 



where near the minima corresponding to the two phases s, b: 
1 



U{n) 



n — ir 



mm I 



(7.4) 



(7.5) 
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The quantities 

TT Z"^ Vn 2 /^, X 

72' T^' 72' ^3, ^4^3 (7-6) 

Us i/s i/s 

are dimensionless functions of x, y and are different in the broken and symmetric phases. 
Including only the quadratic terms the correlator is 

(7r(x)7r(0)) = -^expi-Zv2\^\). (7.7) 
47r|x| 

However, it is more convenient to define a plane-averaged field 

^(^) = /^:J^(^'?/'^) (7-8) 

satisfying 

{niz)nm = J—-^expi-Zv2z). (7.9) 
Area /f 2 

Determining the exponential decay and the magnitude of plane-averaged correlator 
thus gives the parameters Z and V2 of the effective action. 

By scaling vr = (/?///2a)i?^ + const., these formulas can be directly rewritten for R^. 
For a lattice with the geometry N"^ x L the plane Higgs variable is defined as Rp{z) = 
Ea:,{/ -R^(^) ?/) -2)]/^^- At large distances the plane-plane correlator has the asymptotic 
form 

{Rl{z)Rl{0)) = AN-^ie-'""^ + e-(^"-^)™«) + const, (7.10) 

where rriH = Using eq. { \l.9\ ) this determines the quadratic parameters of S'cfr(7r) 

as 

mH = Zv2, Z^ = ]^Al3lmH. (7.11) 

For the potential we write 

U{Tx)a^ = U{^R^ + const.)a^ 

= V2{R' - RlJ' + Vs{R' - RlJ' + V,{R' - Rl^Y + ..., (7.12) 

where the couplings are related to the previous ones by 



2 



2gl\2 

1 x3 



^-Pnyxagl (7.13) 



V3 = V3 \^-Ph 

n = «.ft^(i/3«)' X -L, (7.14) 
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Figure 19: The wave function normalization of the scalar field near the phase transition 
for m*fj = 60 GeV and (3g = 8, 12. The horizontal lines are constant fits to the 
symmetric phase data. The dashed curves labeled v'^/g'^T* connect the MC datapoints 
for that quantity. 



The wave function normalization. The lattice data for the quantity Z'^/g^ is 
shown in fig. ^ for m*^ = 60 GeV and for T near T*. One can see that Z depends 
considerably on T in the broken phase, but is very accurately constant in the symmetric 
phase. The value of Z^ in the symmetric phase is 

= 0.320(5)^73', (7.15) 

where the result is the average of the constant fits to the (3g = ^ and 12 data. In the 
broken phase the behaviour of Zf, can be easily understood via the tree level relation 



The quantity v'^/{glT*) is also plotted in fig. |T9|, calculated from eq. (f4.2|) . We can 
observe that the relation ( [7.16| ) is valid to good accuracy. 

The potential. The potential part U (vr) of the effective theory can be derived from 
the probability distribution p{B?) using the relation ( |7.2| ) and by fitting the parameters 
and Vi in (|7.12|) p|. Using again the m*^ = 60 GeV, /?g = 8 data we are able to 
determine the three parameters -Rmin, ^2 and V3; the statistical accuracy does not allow 



one to fix V4 with reasonable precision. The results of the fits are shown in fig. 2C. 



^This method is known as the constrained effective potential |36 . 
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Figure 20: The lattice results for the effective potential coefficients, determined from 
rriH = 60 GeV, (3g = 8 systems. 



There is an important consistency check for the validity of the effective theory de- 
scription: using eqs. ( |7.1lD and ( [7. 131 ) the mass of the scalar particle is 



m 



H 



8V2Z^ 



Substituting from ( [7.11[ ) this becomes 



4A{^H/a)V2 = l. 



(7.17) 



(7.18) 



Here ^h/cl is the Higgs correlator in lattice units. Eq. ( [T.lSj ) establishes a non-trivial 
connection between the different methods of mass determination. The first method is 
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Figure 21: The lattice results for the test of the effective theory description. 



the direct measurement via correlation functions and the second one is through the 



potential curvature. The lattice data for AA{^h /a,)V2 is presented in fig. It is seen 
that the prediction of eq. (|7.18|) is satisfied reasonably well in both phases. 

Let us turn now to the discussion of the behaviour of the coefficients V2 and V3. When 
the temperature decreases (increases), V2 in the symmetric (broken) phase rapidly de- 
creases, in fact repeating the behaviour of the scalar mass from the correlation function 
measurements. This gives another method for the metastability range determination. 
We write 

oc \T - Ti\'^' (7.19) 

and fit Tj and 7^ separately for symmetric and broken phases. In the symmetric phase, 
the results are 

(3g= S: T_ = 137.0(5) GeV. 7s = 0.51(14) 

/3g = 12 : T_ = 137.3(5) GeV, 7, = 0.47(16). ^ ' 



These values are quite consistent with the numbers in Sec. [7.2| . However, in the broken 
phase we could not obtain acceptable fits. 

It is interesting to note that the coefficient V3 has a strong dependence on temper- 
ature in the symmetric phase. This suggests that the effective potential U{tt) has a 
scale- invariant form in the vicinity of the lower spinoidal decomposition point. If we 
introduce the dimensionless variables 

n = (n - TT^in)/{Zy/m]^), x = xtuh, (7.21) 



36 



2.0 
1.0 
0.0 

K -1.0 

-2.0 
-3.0 



m* = 60GeV (3^ = 8 



► 40 x80 

30^60 



-4.0 
130.0 




135.0 140.0 
rVGeV 



145.0 



40.0 



30.0 - 



20.0 



10.0 



0.0 



m* = 60GeV P^ = i 



-10.0 



MO x80 

•■ 30'x60 



130.0 135.0 140.0 

rVGeV 



145.0 



Figure 22: The temperature dependence of the effective couphng constants and /i4. 



then the action can be written in the form 
S = I d?x 



{dixf + ^T^'' + h-iTl'' + KtI^ + . . . 



(7.22) 



where the dimensionless couphngs hj, and /i4 are expressed via lattice observables as 

/;,3 = V^(2Afl'^, /i4 = 161^2^4^^ (7.23) 
Now, if /is and /i4 do not depend on the temperature, then the system has scale-invariant 



behaviour near the point of absolute met ast ability. The lattice data is shown in fig. |22. 
Indeed, /is is reasonably temperature independent, while the quality of the data does 
not allow to make a definite conclusion on the higher interactions. The values of the 
coupling constants are completely different in distinct phases; they are larger in the 
symmetric phase, showing that interactions are strong there. 



8 Simulations with field, m\ = 80 GeV 

In 0] Monte Carlo studies of the = 80 GeV system, using a lattice action which 
includes the adjoint Higgs field Aq, were reported. The Aq field, a remnant of the 
timelike component of the gauge field, is a heavy field of mass ~ gT, and has been 



integrated out in the theory of eq. ( |1.1| ). The lattice action is 



S = 5[f/,<l>] + i/3GE[Tr^(x)f/;(x)Ao(x + i)[/,(x)-TrA2(x)] 
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Figure 23: The p{B?) distributions for mj^ = 80 GeV, (3g = 12 lattices. The histograms 
are evaluated at I3h = 0.34771. The double-peak structure of the histograms becomes 
less pronounced when the volume becomes larger than 32^. 

+ /32^E^TrA^(x)+/3i^5:[iTrA^(x)]^ (8.1) 
- ^/3HE[^TrA^(x)iTr$t(x)$(x)] 

where 5'[?7, $] is the action given in eq. (|2.4| ), and the lattice field Aq is given in terms of 
the continuum field as Aq = ^g^aTaAQ. The parameters and can be written 
in terms of the 4d continuum variables (see eqs. (11.15-19) and (11.97-99)). 

In lattices of sizes up to 32^, with Pg = 12 and 20, were studied. Clear double- 
peak histograms were observed for both values of signaling a first order transition. 
However, since the publication of 0] we have performed simulations using lattice vol- 
umes up to 48^. Surprisingly, with increasing volume the double-peak feature of the 
histograms becomes less pronounced, in contrast to the expected behaviour in first 
order transitions (see fig. The most striking finite volume effect is the shift of 

the peak positions towards each other when the volume is increased. Qualitatively 
similar finite volume behaviour occurs also in several other systems, which can exhibit 
either first or second order phase transitions in the infinite volume limit (for example, 
the 2-dimensional 7-state and 4-state Potts models, which have respectively first and 
second order transitions ||3y, ^). Thus we cannot yet make definite conclusions about 
the order of the phase transition. 
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Figure 24: The Higgs and W masses of mj^ = 80 GeV, (3g = 12, and V = 40^ lattices 
around the transition temperature. 



However, a note about the computer simulations is in order: the action ( ^.1|) used 
in these simulations includes the Aq field, and because of the Aq hopping term in the 
action it is difficult to write an efficient overrelaxation algorithm. The algorithm we 
use here uses heat bath and Metropolis updates, and it performs more than an order of 
magnitude worse than the overrelaxation described in Sec. ^. Therefore, the statistical 
accuracy of the m*fj = 80 GeV results is not nearly comparable to the m*^^ < 70 GeV 
results, described in Sees. ^|-^. There are also more systematic uncertainties than in 
the simulations in Sees. ^-|^, since the relations of lattice and MS schemes were not 
fully known at the time of the = 80 GeV simulations. 

In fig. ^ we show the Higgs and W masses around the transition temperature. To 
the accuracy of our data, we do not observe any discontinuities at the transition, as 



opposed to the m*j^ = 60 GeV masses in fig. |I6[ However, the errors here are about a 



factor of 5 larger. In the immediate vicinity of the transition, the ratio of the W and 
Higgs masses is my/ Imn ~ 4. 

We measure the critical coupling [3h,c for each individual volume as described in 
Sec. ^. The infinite volume extrapolation is done using volumes larger than 20^, with 
the results [3h,c = 0.347715(8) for (3g = 12, and I3h,c = 0.341700(10) for /So = 20. In 
contrast to simulations with m}^ < 70 GeV, in these simulations was not kept 
constant, but instead the coupling constant Pr was fixed to values Pr = 0.00124 
and Pr = 0.000712, respectively. Using the relation m*^/m^r = ^PrPg/ Ph^ these 
correspond to m*^^ = 79.98 and 79.62 GeV. The transition temperatures and the linear 
extrapolations to the continuum limit are given in table |^. The numbers are somewhat 
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different from tliose presented in f^. Tlie discrepancy is due to an error in tlie 2- loop 
term of the constant pliysics curve in 0. 



/3g 


m^/GeV 


T;/GeV 


12 
20 

oo 


79.976(1) 
79.614(2) 
79.072(6) 


170.0(2) 
169.8(7) 
169.4(17) 



Table 9: Results from ^ 80 GeV simulations, together with the extrapolations to 
the continuum limit. The perturbative critical temperature corresponding to = 
79.072 GeV is about 171.5 GeV. 

The continuum limit T* is somewhat below the perturbative value; however, the sta- 
tistical errors are too large for the difference to be significant. To resolve the situation 
at m*fj = 80 GeV requires simulations performed without the Aq field. Let us note that 
the value of x corresponding to mj^ = 80 GeV is Xc = 0.1188, and the physical Higgs 
mass in the Standard Model corresponding this Xc is about 81 GeV. 



9 What happens at large Higgs masses 

The phase transition in the 3d SU(2)+Higgs system becomes weaker and weaker with 
increasing scalar self-coupling x ~ mj^. An important question is whether the transi- 
tion for some value of x terminates in a 2nd order transition or whether it continues 
as a 1st order one for all x. 

The main difficulty in resolving the order of the phase transition at large scalar self- 
coupling is the presence of the two distinct physical mass scales (inverse correlation 
lengths, m(T) = 1/^). In fact, inspection of fig. shows that for = 60 GeV near 
the transition point, 

^h{T) _ mniT) _1 

/m\ Isymmetric ~ /m\ Ibrokon ~ „• W"*-/ 

mwi-L ) mwi-L ) o 

Furthermore, fig. ^ indicates that this ratio is even smaller for larger m|^, and data 
for smaller indicates that it is larger for them. Note also that vector masses do 
not seem to go to zero at the points of absolute metastability defined in Sec. ^ 

Given the observed hierarchy > ^w, the lattice spacing must be much smaller 
than the vector correlation length, while the lattice size must be much larger than 
the scalar correlation length. This puts stringent constraints on the lattice size and 
makes definite numerical conclusions more and more difficult with increasing rriH (see 
the discussion of the simulations with = 80 GeV in Sec. ^). 

However, if the ratio mniT) /mw{T) continues decreasing with increasing mH-, a 
solution to the problem immediately suggests itself: the vector degrees of freedom 
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become relatively heavy, decouple, and should be integrated out. From the effective 3d 
SU(2)+Higgs theory one thus at large mu proceeds to an effective 3d scalar theory with 
just one degree of freedom. The transition to the scalar theory cannot be carried out 
by integrating out Ai perturbatively, and we first have to discuss the relevant degrees 
of freedom in more detail. 

The 3d Euclidean SU(2)+Higgs theory originating from dimensional reduction corre- 
sponds to a (2+l)-dimensional Minkowskian SU(2)+Higgs theory. This theory has an 
analogy with QCD: we have doublets of scalar "quarks" bound together by an SU(2) 
triplet of "gluons". We have chosen a Minkowskian language in order to keep the 
analogy with QCD transparent. 

Then, one would expect that the lowest lying particle states in the symmetric phase 
are the spin zero scalar bound state and the three degenerate vector bound states, with 
quantum numbers corresponding to the composite operators 

vr = 0V, (9-2) 

= z{(P^D,4> - W; = z{(PW,4^ - {D,4>y4>), Wf = {W+r, (9.3) 

correspondingly. Here = 2X20*. Note that these operators provide also a gauge- 
invariant description of the Higgs particle and intermediate vector bosons in the broken 
phase. In other words, the particle degrees of freedom are the same in both phases. 
This fact alone suggests that there is no distinction between the two phases, so that 
in some region in the parameter space [p^, |38| there may be no phase transition at all. 



However, it does not exclude the situation that going from one phase to another is only 
possible through a phase transition. 

The most general renormalizable local scalar field theory in 3d has the action 

L = ^id,xf + Pix) (9-4) 

where P{x) is a sixth order polynomial which contains five constants (a linear term 
can always be removed by the shift of the scalar field)^ 

^(x) = E-x^ (9.5) 

i=2 ^ 

One of the constants (say, C4) fixes the scale, and the other four (say, C2/C4, Cg/c^, C5/C4 
and Cg) are dimensionless numbers completely fixing the dynamics of the theory. If the 
assumption on the effective theory is correct, then all the four parameters are some 
functions of the dimensionless ratios x and y characterizing the SU(2)-|-Higgs theory. 
In the vicinity of the phase transition, the first of the ratios (C2/C4) is fixed, and there 
is a mapping of the three last ones to the single parameter x of the gauge-Higgs system 

^Note that there is no symmetry (discrete or contimious) which forbids odd powers of % ^ Z^^tt. 
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at T = Tc- Clearly, these relations cannot be computed in perturbation theory, since 
at large the phase transition occurs in the strong coupling regime. After they have 
been found by some nonperturbative means, resolving the order of the phase transition 
on the lattice at large mu becomes a much easier task, because only one essential scale 
is present. 

Finally, note that the action proposed here has a different number of degrees of 
freedom, different symmetries, and a different range of validity than the scalar action 
used in the simulations of [|14|. 



10 What did we learn from the lattice simulations? 

The non-perturbative information we acquired by doing lattice simulations can be 
confronted with different hypotheses on the nature of the phase transition existing in 
literature. It can be used as a test of validity of a number non-perturbative approaches. 



such as the e-expansion the renormalization group approach and the Dyson- 
Schwinger equations [|^. In this section we provide a summary of our findings and 
compare them with some of the mentioned approaches. 

10.1 Comparison with perturbation theory 

Perturbation theory can be used for the computation of different quantities in the 
broken phase. Not many quantities are known to high order in perturbation theory 
in the 3d SU(2)+Higgs theory, though. The whole list consists of the 2-loop effective 
potential for the scalar field in different gauges p|, ^ (and the values of different 
condensates related to the effective potential, see 0), and the 1-loop pole masses of 
the gauge and Higgs particles (see |H2| and below). 



The effective potential. If the temperature is fixed to some value, the vev of the 
Higgs field, defined by eq. ( ^.4|) , can be computed with the use of the 2-loop effective 
potential. At the same time, the "exact" vev of the Higgs field can be found on the 
lattice (after appropriate extrapolation to the continuum limit). The 2-loop effective 
potential gives a reasonably good prediction for v/T. For example, 5v/v ~ 0.3% for 
the vev of the Higgs field when v/T ~ 2, and 5v/v ~ 10% for v/T ~ 0.64 (these 
numbers refer to m*^ = 80 GeV). Here 6v = t'2-ioop — "^exact- The exact vev is found 
to be smaller than that from 2-loop perturbation theory. This proves that the 3-loop 
linear correction to the effective potential comes with a positive sign. If the term in 



eq. ( ^I3|) is added to the 2-loop effective potential with /5 = 50, then the agreement 
between the exact and 3-loop results is within 1% even for v/T as small as 0.6. In 
other words, the temperature dependence of the vev of the Higgs field can now be 
found analytically with 1% accuracy up to rather small values of v. This result is even 
more non-trivial than a direct analytical 3-loop computation of the effective potential, 
because it indicates that higher loop corrections are indeed small in this region of vevs. 
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The determination of 3-loop corrections to the effective potential allows one to get 
a better idea of the expansion parameter of the SU(2)+Higgs theory. As we have 
discussed in it is difficult to get a reliable estimate of the expansion parameter for 
the effective potential on the basis of 2-loop perturbative computations. The expansion 
parameter (see Sec. |10.1|) is proportional to gl/^nniT) = ^g^./ijcj)), but how big is the 
constant of proportionality? Indeed, the loop expansion may be reliable in the broken 
phase in spite of the fact that the 1-loop correction is comparable with the tree term, 
just because the scalar self-coupling constant is small. Moreover, the magnitude of the 
2-loop correction may be changed by adjusting the scale /x, so that the 2-loop correction 
may even be tuned to zero by some choice of /i. Hence the determination of the 3-loop 
correction to the effective potential is important. 

Let us take for definiteness m*^ = 82 GeV, i.e., Xs/g^ = 1/8 according to eq. ( p.9|) . 
In addition, consider the region near (take mg = 0) and include only gauge field 
loops (neglect mi and ^2). Then the estimates of the different contributions to the 
effective potential are 
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where mT((f)) = gz<t>/'^- Inspection of these relations immediately shows that a rea- 



sonable estimate of the expansion parameter is — This quantity is rather large 



9l 

■iTmT{4>) ' 

(~ 0.8) in the vicinity of the phase transition for m*^ = 70 GeV, and the convergence 
of perturbation theory is expected to be quite bad. The main reason why the vev of 
the Higgs field can be found with a pretty high accuracy is that the effective potential 
is known up to a high order in perturbation theory. 

Perturbative correlators in the broken phase. The next two quantities char- 
acterizing the broken phase are the vector and scalar correlation lengths. 

As in Sec. ^, at a finite order in perturbation theory there are different methods of 
calculating the correlators. In principle, the most straightforward is a strict pertur- 
bative calculation in powers of h in the loop expansion. Such a calculation can be 
done by shifting the Higgs field to the classical broken minimum where the Goldstone 
boson mass m2 vanishes, and by then calculating all the 1-loop diagrams, including 
the reducible tadpoles (the calculation in ||2| is organized in this way). This method 
gives an explicitly gauge-independent results for the pole masses. Unfortunately, such 
a calculation gets increasingly unreliable as one approaches the phase transition, since 
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the tree-level results mi, for the Higgs and W correlators, respectively, go to zero. 
The same problem was met in Sec. ^ in connection with the effective potential. One 
should therefore solve the relevant equations numerically, hoping that the higher-order 
corrections included are the essential ones. With this method, a small residual gauge- 
and /i-dependence remains, and its magnitude can be used as an estimate of the con- 
sistency of the approximation. We find that the gauge- dependence of the final results 
for the correlators is numerically small, at most of the order O.OSg^. Below we give 
results in the Landau gauge. The /^-dependence is also briefly discussed. 

To calculate the correlators, we first solve for the radiatively corrected location of 
the broken minimum (pb from the RG-improved 2-loop effective potential V{(j)) 0. The 
requirement of RG-improvement used here for the renormalization scale is that 
the 2-loop contribution to the derivative of the effective potential, d^V2{4>, /u), vanishes. 
With (pf, fixed, we calculate the inverse propagator, and solve numerically for the real 
part of the location of the pole. Since (pj, is an approximation to the exact minimum of 
V{(f)), no reducible tadpole diagrams are needed in the calculation of the propagator in 
contrast to the /^-calculation. The renormalization scale used in the calculation of the 
propagators is chosen to be fi{4>b)- The dependence on /i is formally of higher order 
than the accuracy of the 1-loop calculation, but is numerically visible, see below. The 
details of the calculation are in [Appendix C . 

The pole mass obtained from the correlator for the W particle is shown in fig. [1^. 
In this figure we also show the tree-level value rriT- It is seen that in the broken phase, 
1-loop perturbation theory gives excellent results. This is quite unexpected since the 
expansion parameter is rather large at the transition point. In fig. [16], we show also the 
pole mass of the Higgs particle, together with the tree-level value. There is also a third 
curve: from the known 2-loop contribution V2(0) to the effective potential, one can 
derive the momentum-independent part of the 2-loop self-energy: 11^(0) = — Vg^^b)- 
When this is added to the complete 1-loop result, one gets the dotted curve. It is 
seen that the discrepancy between lattice and perturbation theory is larger than for 
the W correlator. The effect of V2 suggests that a complete 2-loop calculation of the 
Higgs self-energy might give a better estimate. The importance of 2-loop corrections 
for the Higgs particle is in accordance with experience from the effective potential, 
see [0. Finally, let us point out that the tree- and 1-loop W masses, as well as the the 
Higgs mass with the 2-loop contribution included, are practically independent of /i. 
For the tree- and 1-loop Higgs masses, the /i-dependence in varying /i in the range 
O.5/i(0b) . . . 2.O/i(0fe) is of the order of 0.03c/|. 

Attempts to describe the correlation lengths in the symmetric phase can be found 

Parameters of the phase transition The characteristics of the phase transition, 
such as the critical temperature, latent heat, bubble nucleation rate, surface tension, 
and correlation lengths in the symmetric phase cannot be defined in the perturbative 
framework only, because the symmetric phase is in the strong coupling regime. The 
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failure of perturbation theory is clearly seen in order by order computations. For exam- 
ple, the vector 4-loop contribution has a logarithmic singularity at = 0. Nevertheless, 
it is still interesting to compare 2-loop and 3-loop predictions with the results of lattice 
simulations. The RG-improved 2-loop effective potential is real in the vicinity of the 
critical temperature for the whole range of and is regular at the origin. Hence, the 
computation of the above-mentioned characteristics of the phase transition is not faced 
with formal mathematical difficulties. 

The critical temperature of the phase transition derived by 2-loop perturbation the- 
ory is presented in table |^. It is somewhat larger than that derived on the lattice. 
The vev of the scalar field at T* for m*j^ = 60, 70 GeV coincides with the lattice value 
within errorbars. This perfect agreement, however, is an incident, since the comparison 
of vevs is done at different temperatures. The same is true also for the latent heat, 
since it is mainly determined by the value of the scalar condensate in the broken phase. 
We do not expect the same coincidence at larger Higgs masses. 

One of the crucial quantities for the computation of the bubble nucleation rate is 
the surface tension. For m}j = 60 GeV, its value on the lattice is considerably (about 
3 times) smaller than that derived from perturbation theory as 



Note that this equation is valid to leading order only; in higher orders the wave function 
renormalization must be taken into account. However, the huge discrepancy with the 
lattice result shows that perturbation theory is not applicable at all for the computation 
of the surface tension, at least for m|^>60 GeV. 

To summarize, the computation of the characteristics of the phase transition from 

2- loop perturbation theory (for 35 GeV< <70 GeV) is accurate for the critical 
temperature within 0.8 — 1.6% and for the ratio v/T within 6%. Perturbation theory 
fails to describe the surface tension (and, therefore, the bubble nucleation rate) at least 
for m*jj = 60 GeV. 

Adding to the effective potential the 3-loop linear term determined in Sec. ^ does 
not improve the accuracy of perturbative predictions. We take as an example m*^ = 60 
GeV and /3 = 50, as follows from lattice simulations. The critical temperature remains 
roughly as far from the lattice T* as at 2-loop level, being now smaller than the lattice 
value, T3-i°°p = 136.8 GeV. The 3-loop ratio ^;3-ioopj2^3-ioop)/2.3-ioop ^^^^^ 20% larger 
than the exact value[3. The 3-loop perturbative latent heat is away from the lattice 
results by almost 50%. The most drastic deviation is in the surface tension. The 

3- loop value of it is a/{T*Y — 0.018 - a factor 7 larger than the lattice number! 
The behaviour of the effective potential in two different situations is illustrated in 

^"^The estimate of v/T in the broken phase at a given temperature of course does improve, since the 
3-loop part is determined by just this requirement. The discrepancy is due to the fact that the exact 
and 3-loop critical temperatures are different, and v/T is strongly temperature dependent near Tc- 
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Figure 25: The RG-improved 2-loop effective potential (dashed line) and the 3-loop 
effective potential with /5 = 50 (sohd hne) for = 60 GeV at the non-perturbative 
critical temperature T* = 138.38 GeV. Note that in this figure we have expressed in 
4d units [04d = 03d(T*)i/2]. 

figs. |2^, |2^. It is clearly seen that perturbation theory cannot quantitatively describe 
the phase transition. 



10.2 Comparison with other non-perturbative approaches 

We were not able to make a detailed comparison of the results of lattice simulations 
with all the non-perturbative approaches suggested in literature. For example, in p9[ 
devoted to the application of the e-expansion, the authors compare their results to 1- 
loop effective potential predictions. As we argued in p[, 1-loop computations in 3d have 
a considerable unphysical scale dependence, and this fact makes a direct comparison 
difficult. We could not make a comparison with the renormalization group approach 

I, m 



m 



40| , either, nor with |^. Hence we discuss here |^ and 



In 1^ it was suggested that the non-perturbative effects may considerably modify 
the effective potential of the Higgs field near the origin. Assuming that perturbation 
theory works in the broken phase, a measure of the non-perturbative energy shift at 
= was introduced, 

12 



9l 
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Figure 26: The 2-loop and 3-loop effective potentials at the corresponding perturbative 
critical temperatures, T*^ = 140.25 GeV and 136.82 GeV, respectively. 



Here T* is the exact value of the critical temperature. It was conjectured in |^ that 
if Ap is positive then the phase transition is stronger than predicted by perturbation 
theory. In particular, if is sufficiently large (say, Ap = 0.4), then the lower metasta- 
bility temperature is considerably smaller than the critical temperature, and the vev 
of the Higgs field at T_ is substantially larger than it is at the critical temperature. 
The lattice simulations allow to check the validity of this hypothesis. Let us take for 
definiteness = 60 GeV. Then, Ap is to be found from 3-loop perturbation theory 
(2-loop perturbation theory gives a result off by 20% for the vev of the Higgs field and 
therefore is not to be used). We get Ap ~ —0.08, see fig. |25|^^| . In other words, the 
phase transition is weaker than predicted by 3-loop perturbation theory. Hence the 
assumption of the dominance of the gluonic condensate contribution to the vacuum 



energy in the symmetric phase, used in [^], appears not to be satisfied . 

In P2| it was suggested that a reasonable description of the symmetric phase can be 
achieved with the 1-loop Schwinger-Dyson equation. From the analysis of this equation 
it was anticipated that the symmetric phase can be interpreted as a Higgs phase whose 
parameters are determined non-perturbatively. One of the predictions of this approach 
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'^In we had another estimate of the parameter Ap, Ap = 0.2. The difference is because in [|| 
2-loop perturbation theory was used, the 2-loop relations of the lattice and MS schemes were not fully 
available, and the proper extrapolation of the lattice results to the continuum limit was not done due 
to the limited amount of data. 
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is that the vector boson mass in the symmetric phase is smaller than 0.27gl, what is 
to be confronted with the lattice value mw ~ l-lfi'i- We believe that the main reason 
why this approach does not work is because it relies on an accurate description of the 
symmetric phase by 1-loop perturbation theory. In particular, it was estimated in ^ 



2 

that the expansion parameter in the 3d theory is tt--^, while we find it to be about 6 
times larger. 

In ^ it was suggested that the introduction of a "magnetic mass" rriM ~ S'l/Svr 



in the propagator of the gauge boson will cure the infrared problem and allow one 
to estimate the magnitude of non-perturbative effects. This recipe would effectively 
produce a linear 3-loop term with a negative coefficient: 

~ -^{mli + m^Tf = + IgsMli + •••)' (10-4) 

where ttit = gz<p/2. This is clearly in contradiction to the existence of the positive 
linear term we found in the simulations. 



11 Results for different physical 4d theories 

So far we have given results for the characteristics of the phase transition in the 
SU(2)+Higgs model using the simplified relations of 4d and 3d parameters given in 
eqs. (|2.8| )- (|2.1CI|) . Here we explain how these results can be converted into character- 



istics of finite-temperature phase transitions in physical theories where the relations 
between 4d and 3d are more complicated. We also give explicit results for the phase 
transition in the 4d SU(2)-|-IIiggs theory without fermions, in the 4d SU(2)-|-Higgs 
theory with the fermionic content of the Standard Model, and in the Standard Model. 

To begin with, let us restate the observables relevant for the phase transition in 
terms of properties of the 3d theory alone. The parameters of the 3d theory are 
gl,x,y, defined in eqs. (|2.8|) - (|2.1CI|) . From lattice measurements, one can derive the 
expectation value of dimensionless gauge-invariant observables, like {(p'^ 4>{gi)) / gl (one 
must use eq. { ^^) to change scheme from lattice to MS at = gl). Let us denote this 



particular expectation value by £3, £3 = {(j)^ (pi^gD) / g^. Let l\ be the value of £3 in the 
broken phase at the critical point (xc, i/c), and A£3 be the difference of the values of £3 
in the broken and symmetric phases at the critical point, Al^ = ^3 ~ ^3- The critical 



point {xc, Uc) is determined as explained in Sec. pA\ The observable Ai^, is related to 
the dimensionless quantity e3{x,y) defined by 

exp[-V3gle3{x,y)] = J V4>VAexp{-S[cP, A]), (11.1) 

through 

^^IpllA = Ms (11.2) 

dye 
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where Ae3(x, = el{x,y) — el{x,y) and V3 is the volume of the system. Finally, we 
define the dimensionless quantity (T3 related to the surface tension by 



1 p^^^ 

(73= hm — r-rln^T— , (11-3) 



in analogy with eq. (pl5|). The observables Xc, Vc, ^3, ^^3 and are all dimensionless 
quantities defined strictly inside the 3d theory, are independent of the parametrization 
used in eqs. (|2.8|)-( prT0D , and are measurable, in principle, to arbitrary accuracy with 
lattice simulations, since the lattice counterterms of the 3d theory are known exactly. 
Next, one needs the relations to 4d. The relations to 4d consist of two parts: 

(a) One needs the expressions of the parameters g^, x and y of the 3d theory in terms 
of temperature and the physical parameters of the 4d theory. For definiteness, let us 
take as one of the physical parameters the pole mass of the Higgs field, mn- Then 
the parameters of the 3d theory are of the form g^ = g^{mH,T), x = x{mH,T) and 
y = y{mH,T). 

(b) One needs the relation of 3d observables to 4d thermodynamics. This relation is 
completely determined by the equation p{T) = —Tg^e^{x,y), which holds apart from 
inessential mass-independent terms. 

Using the relations in (a) and (b), one can relate Xc, ?/c, ^3, ^£3 and to 4d quantities 
characterizing the thermodynamics of the phase transition. First, from the equations 

I = ^/^^'I^^,) (11.4) 

one can solve for the Higgs mass rriH and critical temperature Tc to which the phase 
transition in the 3d theory at the point {xc, yc) corresponds. When Tc is known, one 
can calculate the gauge coupling: g^ = g^^mn, Tc). Then one can determine v^{Tc)/T^, 
defined by v'^{Tc)/T^ = 2{(j)\(j)^{Tc)) /Tc according to eq. (p. 4]), from 



T2 Tr 
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' 1671^ Tc 
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The surface tension cr/T^ is obtained from 



^ = ^3^. (11.6) 



For the latent heat L = Tdp'^iT) — p[(T)], one gets 

L_ _ J__d_ 
T^ ~ T^dT 



TglAe3ix,y) 



T=Tc 



At4(%-^'^]. (117) 
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Here we utilized eq. (|11.2|) together with the fact that at the critical line, the following 
two equations hold: 



Ae3(xe,i/c) = 0, (11.8 
dAe3{x,y) dycdAe3{x,y) 



dx dxc dy 



11.9) 



Note that both dx/dT and dy^dxc in eq. ( |1 1 . T] ) are non-zero through loop corrections: 
dx/dT is no n- vanishing due to logarithmic 1-loop corrections to the dimensional reduc- 
tion of the coupling constants (73, A3, and dy^dxc is non-zero due to loop corrections 
inside the 3d theory (at tree-level, the phase transition takes place at ?/ = independent 
of x). Numerically, dyddxc ~ —1 but dx/dT <^ dy/dT, so that in realistic cases 

{dx/dT){dyJdxc) m in^ 
d^ ^^^-^^^ 

Hence one does not need to determine the derivative of the critical line y = ydx) with 
as good a relative accuracy as the jump Ai^ of the order parameter at the point (xc, yc)- 
If one uses the parametrization of eqs. (|2.8| )- (p.lCI|) in eq. (ITLTD , then dx/dT* = and 
one gets the expression in eq. (|5.3|). 



To get the values of £3, Ai^ and o"3 from the results given in Sec. ^, one has to use 
eqs. ( [11.5| )- ([11.71 ) in the inverse direction, employing the parametrization of eqs. (|2.8| )- 
( p.lOj ). For instance, the surface tension cr /T*^ = 0.0023 for = 60 GeV given 
in table corresponds to 0-3 = 0.0023/0.44015^ = 0.0119 according to eqs. (|1L|) 
and (|2.8| ). The values of Xc, yc, ^3, Ai^ and obtained this way are shown in the 
second block in table |TU[ To go back to 4d units for different physical theories, one 
needs the true values of g^/Tc, x'{Tc) and y'{Tc). In general, the value of g1/Tc differs 
from 0.44015, so that cr /T'^ differs from the value in table ^ see table |10. 

The explicit form of the expressions for g^, x and y in terms of temperature and the 
physical 4d parameters of the Standard Model to order g^ have been given in , using 
the approximation g''^ ~ g'^. With these relations, we can give results for the thermo- 
dynamical properties of the EW phase transition. To be more precise, we will consider 
three different theories reminiscent of the EW sector of the Standard Model. As a start- 
ing point, we consider the case mz = my/ so that g' = 0. This is a SU(2)-|-Higgs theory 
with the fermionic content of the Standard Model. The functions x{mH,T), y{mHiT) 
and gl{mH, T) in this theory, with mtop = 175 GeV, have been given in figs. 7-8 of 0. 
The results for the phase transition, obtained from the values in the second block in ta- 
ble |10| using eqs. ( |11.4| )-( p!l.7| ), are shown in the fourth block in table |10[ Results for the 



full Standard Model can be obtained from the SU(2)-|-Higgs-|-fermions model by taking 
into account the U(l)-subgroup perturbatively with the help of fig. 9 in The results 
are in the bottom block in table |10[ Finally, we give results also for the SU(2)-|-Higgs 
model without fermions. This case is obtained from the SU(2)-|-Higgs-|-fermions model 
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by putting gy = np = as = and fixing mw = 80.22 GeV. The results are shown in 
the third block in table 0. In principle, the results for the SU(2)+Higgs theory should 
allow a comparison with the 4d lattice simulations in |2^, At present such a 
comparison is not rigorously possible, however, since the relation of the gauge coupling 
gj^ used in 4d simulations to continuum physics is not known. 



12 Applications to cosmology 
12.1 The phase transition 

The results of non-perturbative lattice MC simulations allow one to considerably re- 
duce the uncertainties associated with the dynamics of the electroweak phase transition. 
Strictly speaking, they still refer to a somewhat unphysical situation, since the U(l) 
subgroup present in the electroweak theory is omitted. As we discussed in the previous 
section, the perturbative effects associated with the U(l) subgroup are rather small; 
the estimate of the non-perturbative effects requires lattice simulations in the com- 
plete 3d SU(2)xU(l)+Higgs model. All specific estimates in this section are based on 
the assumption (quite reasonable, though) that all perturbative and non-perturbative 
effects of the U(l) subgroup are small. 

Let us reconstruct the picture of the phase transition for = 60 GeV, for which 
we have the best lattice data. For the MSM with the top mass rrit = 175 GeV this 
corresponds to the pole Higgs mass of mn = 51.2 GeV, excluded experimentally, but 
for a number of extensions of the Standard Model this may be quite realistic. To make 
the discussion less model dependent we will use the variable T*; with the help of the 
results of and Sec. |TT| everything can be re-computed for any specific model. We 
shall omit here the errorbars from the numbers. 

The critical temperature of the phase transition is T* = 138.4 GeV, and the vev- 
temperature ratio at T* is v/T* = 0.67. The scalar correlation length in the symmetric 
phase is ~ and in the broken phase C,b ~ The corresponding vector 

correlation lengths are a factor of three shorter. The domain wall separating the 
broken and symmetric phases has the surface tension a ~ 0.002T^*'^, and the profile of 
the scalar field has asymmetric tails on the different sides of the domain wall, 6(j){x) ~ 
exp(— At x* y. = 138.7 GeV only the symmetric phase is stable, and 
at T* < T* = 137 A only the broken phase is stable, while in between both phases 
can exist simult aneouslyp^. The bubble nucleation temperature Tbubbie lies somewhere 
between Tc and T_ and may be estimated with the use of the surface tension and latent 
heat found on the lattice. 



^^We use here the results of subsection 7.2 which are the most accurate 
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Table 10: Properties of the phase transition in different physical 4d theories. The 
simulations with = 35 GeV do not correspond to any physical pole Higgs mass 
in the SU(2)+Higgs+fermions model, at least according to the 1-loop formulas for the 
couplings used in f^. The errors for the 2-loop perturbative results in the first block 
indicate the effect of varying in the range 0.5/iopt • • • 2.0/iopt- The error estimates 
for the 3d lattice SU(2)+Higgs model follow from the error estimates in Sees. |5.1[ - |5.3| . 
The surface tension measurements for = 35 GeV are in parentheses for reasons 
explained in Sec. |5.3| . In the 4d results, we only show the central value. For the 4d 
SU(2)+Higgs model, the additional relative error from the relations to 4d should be 
below 1%, and for the SU(2)+Higgs+fermions model, it may be a few percent (in the 
critical temperature, the error is an order of magnitude smaller) f^. The last block 
concerning the Standard Model is based on a purely perturbative estimate of the effect 
of the U(l)-subgroup fig. 9], and the errors may be large if the non-perturbative 
effects related to the U(l)-sector are significant. 
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The thin wall bubble nucleation rate is given by 



r = «:T^xp(-^^) (12.1) 

where Rc = 2a/ e is the radius of the critical bubble, e = is the pressure difference 
between the broken and symmetric phases, AT = Tc — Tbubbie, L is the latent heat 
of the transition, and k is the prefactor. Estimates for k in different models and 
approximations can be found in |^S|, The bubble nucleation temperature is roughly 
determined from the condition 

~ 1. (12.2) 

c 

Inserting the lattice numbers to these relations gives an estimate ~ 0.001, i.e., the 
bubble nucleation temperature is very close to the critical one, T^^^bbic = 138.3 GeV. 
The smallness of AT*/T* is due to fact that the ratio a^/L'^T* ~ 10~^ is so small ^|. 



Since AT*/(Tc — T_) ^ 0.1 is also small, one is in the thin-wall regime; indeed, the size 
of the bubbles when they nucleate is at least Rc — 110/T^* which is much larger than 
the scalar correlation lengths in the broken and symmetric phases at Tbubbie- Since 
^bubble is very close to the critical temperature, the vev of the Higgs field at Tbubbie is 
almost the same as at T*. 



12.2 The out of equilibrium condition for electroweak baryo- 
genesis 

One of the motivations for the study of the electroweak phase transition is its ap- 
plication to electroweak baryogenesis. The rate of the anomalous baryon number 
non-conserving processes is high in the symmetric phase but is suppressed by 
the Boltzmann exponent in the broken phase |^|. Baryogenesis occurs at the 1st or- 
der electroweak phase transition, and the mechanism-independent constraint on the 
strength of the phase transition is that the rate of fermion number non-conservation 
in the broken phase at the bubble nucleation temperature be smaller than the rate of 
universe expansion 
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We parametrize the rate of sphaleron transitions in the broken phase as 



2Esph{T) 
^T 



exp 



-^sph(^) 

T 



:i2.3) 



where the factors A^tr — 26 and A^rot — 5.3 x 10^ are zero mode normalizations |Q and 
Esph(T) is the effective sphaleron mass at temperature T. Then the out of equilibrium 
constraint reads 

-E'sph(rbubble) /^bubble > 45. (12.4) 
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In principle, this bound can be converted into an upper bound on the parameter Xs/g^, 
completely defining the dynamics of the 3d theory. An exact determination of this 
bound would require a non-perturbative evaluation of the sphaleron rate in the broken 
phase (for lattice simulations of topology changing processes in the symmetric phase, 
see 0). Unfortunately, this problem has not been solved yet. Below we will esti- 
mate the critical value of Xs/g^ assuming that the expansion parameter in the broken 
phase is gl/{7imT) = 2g3/{n(f)b) - a value inspired by the higher order computations 
and lattice simulations. 

The best estimate of the sphaleron rate available now is the 1-loop computation 
in where the determinant of small fiuctuations was computed numerically in a 
bosonic theory]^. According to |5^, the rate in 1-loop approximation is just given by 



eq. (|12.3| ), where 



T gi g^ T 

and 0b is the vev to be determined from the 1-loop effective potential (measured here 
in 4d units). Arguments in favour of absorption of the 1-loop effects into the vev of the 
scalar field have been presented in |5^. Some numerical values of the function B are: 



for = 0, 40, 45, 50, B = 3.04, 3.41, 3.44, 3.48, correspondingly [g. We assume then 
that the exact value of the sphaleron mass is given by ( |12.5| ) plus corrections, and that 
(pb/T is to be replaced by the exact gauge-invariant value v/T determined on lattice. 
Then the first correction is of the order A{g'^/{7TmT))'^, where A is a number of the 
order of unity. A conservative limit on the ratio v/T is obtained when A is positive, so 
that v/T > 1.22 (we take B corresponding to = 50 GeV). At this value the 2- loop 
correction is about 10% and the 3-loop correction is expected to be of the order of 4% 
and can be neglected. If, in the contrary, A is negative, we get v/T > 1.49. 

Now, we choose the weaker constraint v/T > 1.22 and convert it into an upper limit 
on m*fj. Since we do not have lattice simulations for the whole range of values 
from 35 GeV to 60 GeV, we take the 2-loop predictions for v/T; as we discussed, this 
is accurate within a few percent. We find that if m*jj ~ 42 GeV then v/T^ = 1.22. The 
bubble nucleation temperature is somewhat smaller than the critical temperature. The 
thin wall approximation for the tunneling rate is not applicable here. Assuming that 
the perturbative description of bubble nucleation is valid in this region of the Higgs 
masses, we estimate that v/T aX the nucleation temperature is about 20% larger than 
at the critical temperature. This is derived by defining the bounce for the action 

^(50)' + V^2-loop(0) (12.6) 



^■^ Recently the fermionic determinant in the background of a sphaleron was computed in pT\ . 
The authors concluded that the fermionic contribution suppresses the rate and is numerically very 
important. However, the effect of fermions can be absorbed into the definition of the 3d coupling 
constants; after this the fermionic contribution is negligible, see |5q|. 



54 



numerically, and requiring that the value of the bounce action is ~ 140 (see, e.g. 
For example, for m*H = 43 GeV, T* = 108.6 GeV, v{T*)/T* = 1.16, TbubWe = 106.3 
GeV, w(Tbubbie) /^bubble = 1.43. Then the requirement v/T > 1.22 gives mjj < 46 GeV. 
The perturbative account of the U(l) factor makes the phase transition stronger first 
order, correcting this bound by a factor 



J^(2 + ^^) (12.7) 
V 3 cos"* Uw 

following from the 1-loop effective potential. We then finally get m*jj<50 GeV. If 
the correction to the sphaleron rate is in fact negative, then the number is smaller, 
about 45 GeV. It is interesting to note that this bound is very close to the initial 
1-loop computation in [^. We stress, however, that these estimates are subject to 
verification by the future lattice simulations of the theory with the U(l) subgroup and 
to non-perturbative evaluation of the sphaleron rate. 

From eq. ( |2.9| ), < 50 GeV corresponds to x < 0.043 and m*jj < 45 GeV to x < 
0.034. To summarize, the upper limit to the parameter Xs/g^ in the 3d SU(2)xU(l) + 
Higgs theory is likely to be 

Xs/gl < 0.04. (12.8) 

In order to define the constraints following from this requirement on the particle 
spectrum of the underlying 4d theory, one has to express this ratio through the physical 
parameters of the 4d theory at the critical temperature. This computation may be quite 
involved but it is very clean from the physics point of view and does not contain any 
infrared divergencies. An essential point is that only 1-loop graphs need be computed. 
Indeed, a 1-loop computation provides O(a^) accuracy in the coupling constants of 
the effective theory. Moreover, the critical temperature enters to the ratio ( |12.8| ) only 



through logarithms, so that even a 1-loop estimate of it will give sufficient accuracy. 



The application of the constraint of eq. ( |12.8D to the case of the Minimal Standard 



Model follows immediately from fig. 8 of 0. Indeed, if rrit = 175 GeV, then no Higgs 
mass can ensure the necessary requirement of eq. ( |12.8|) P^. If the top quark were lighter, 
then some low Higgs mass value might be possible, see fig. 



According to |£2[ the phase transition in the MSSM occurs in the same way as it does 
in the MSM. If true, then MSSM also fails in generating a sufficiently strong first order 
phase transition. The two Higgs doublet model has more freedom, and the results of 



indicate that the constraint ( 12. 8| ) can be satisfied there 



^''To be more precise, the computations in relating the physical masses of the W boson and 
Higgs particle to the parameters of MS scheme break down if the physical Higgs mass is close to the 
Coleman- Weinberg limit. In this limit the higher order Yukawa corrections start to be important in 
the procedure of dimensional reduction as well. So, it is not excluded that Higgs masses close to the 
Coleman- Weinberg limit are still possible. 
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Figure 27: The critical value Xc = ^3/9^ as a function of the physical Higgs mass 
niH and the top quark mass nitop. In general, x depends on the Higgs mass, the top 
mass and logarithmically on the temperature. In calculating Xc, we have solved the 
temperature from the equation y = 0, which according to table |10| is very close to 



the true critical value of j/c. The small error in Tc does not affect Xc much due to the 
logarithmic dependence. The value of the dimensionless U(l)-couphng g'^/gl is not 
shown in this figure; it is g'i / g1 ~ 0.3. 



13 Conclusions 



The 3d formalism, developed in the series of papers [jTl-040 provides a powerful 
tool for the study of phase transitions in weakly coupled gauge theories. It factorizes 
the perturbative and non-perturbative physics and allows one to construct effective 3d 
theories, describing in a universal way phase transitions in a large class of underlying 
4d theories. The effective theories in 3d contain bosons only, and may be used for high 
precision lattice Monte Carlo simulations. 

In this paper we reported on lattice simulations in the 3d SU(2)+Higgs model, which 
is an effective theory for the SU(2) sector of the MSM and its extensions. The nature 
of the phase transition at moderate Higgs masses mu < m\Y is clarified, and the 
results presented here form an "experimental" basis for different theoretical schemes 
attempting to describe strong coupling phenomena at T ~ Tc. 

From the phenomenological point of view, to our mind, the most interesting further 
problems to be solved are the role of the U(l) factor in the phase transition, and the 
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rate of the sphaleron transitions in the broken phase near Tc- The solution of the first 
problem on the lattice is a straightforward generalization of the SU(2) case considered 
in this paper. The non-perturbative estimate of the sphaleron rate is a much more 
complicated problem, and it is even not clear how it can be solved in principle. 
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Appendix A 

In this Appendix we rederive the result (11.76-78) for the h (loop) expansion 

in the broken phase of the 3d SU(2)+Higgs theory writing it in a more explicit and 
complete form, and discuss the parametrisation of the 3-loop terms (the statistical 
uncertainties do not allow a determination of the 4- loop terms). The computation 
takes place in three steps: first find the potential (gauge and /i dependent), then find 
the value V{y) of the potential at its broken minimum (gauge independent but /i 
dependent )|3, and finally calculate the condensate from 

The starting point is the potential up to 3 loops {fi is a loop counting parameter): 

V = Vo + Wi + h^V2 + h^Vs. (A.2) 



Perturbative computation gives V in the form V = V{mT, mi,m2, 4>; fJ.), where 

rriT = ^934>, ml = mKfi) + SX^cj)'^, ml = ml{^) + Ag^^. (A.3) 

Thus, equivalently V = V {niKfi) , (j); ft) . Here, in contrast to eq. (|A.1|) , denotes the 
real Higgs field. 

Note that we have simphfied notation here: by v we mean the the location of the minimum of 
the effective potential, previously denoted by 4>b- The v used here should not be confused with the 
v{T) used previously and defined as v^{T)/T^ = 2(0t0(r))/T. 
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The tree and 1-loop potentials are 

V, = + 1a30^ (A.4) 

Vi = -^{6ml + ml + 3ml), (A.5) 
and the 2-loop part in the Landau gauge is in eq. (1.33), in the general covariant gauge 



in and in the -R(C) gauge in . The 3- loop potential is near the classical broken 
minimum of the general form 

^3 = Tr-2T7-^ log — +o -TTT^ 2 + -/1 + ^3^^, (A.6) 



9mil67r2V ruT 2/ 128 V 2 / (47r)3 (47r 
where 

— ^ = - — (mi + 3m2), h = — , (A.7) 
amg 87r 

and /3 is a complicated function of the masses, which we only need near the broken 
minimum. Note that h in eq. (|A.7|) differs from the definition in eq. ( |2.11 ). 



The n dependent terms of V3 are known (eq. (1.70)) since they should cancel the 
/i dependence of Vi so as to make the /i dependence of Vi + V3 to be of order . 

The terms proportional to powers of l/m2 are singular at the tree minimum of the 
potential. They identically cancel when the value of the potential at the minimum is 
computed perturbatively, see below. 

The effective potential itself is gauge dependent but its value in the broken minimum 
is gauge independent (but dependent on the scale n). The minimum is defined by 

ay(m^(/i),0;/i) ^ ^ 

d(j) 

Denoting the solution of this equation by = f one can solve for it in the loop {%) 
expansion: 

- ^^(0) ^ '-(1) ^ - 

with the result 

4) = 4K'Ko))K(^(o)) - f V^2(4)) (A.12) 

3 3 

3 + T^ + — 1 + -^^(^(0)) 



= vl. + Tivl. + h\l. + OiK"), (A.9) 



2 



167r2 



3,2/ 8 \ / 4 \ 5/i ,2 /i^ 3/i^ 



16 V /iVV /iV 4 2 8 2 + /i 
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Here V^^,(f^o)) i^san writing Vn = l^(m|,0), taking partial derivative with respect 



to 0^ and evaluating the resuh at the saddle point value ([A.10| ). Hereby one will need 
the saddle point values of the masses: 



mi = hmx = \J —2m\{ix), m2 = 0. 
The saddle point values satisfy the equations 

drfiT 1 drfii 1 



dm\ 



dm\ 



mi 



(A.13) 



(A.14) 



The pole terms l/m2 in (|A. 12|) cancel before taking the limit m2 0. These arise 



when taking derivatives of Vi ~ m\, V2 ~ "^25^3 ~ l/'^2) etc. In gauge invariant 
quantities calculated below they cancel and will not be explicitly written down. 
Inserting (|A.9| ) to ( |A.2|) and expanding one gets 



V{v) 



V, + nvi + n\v2 - X',\v[f] + k'W^ + {v^ + -v^vl^) 



(A.15) 



where all quantities on the RHS should be evaluated at the saddle point values ( |A.13|) . 
Note that f ^ is only needed up to 1 loop. 

The potential at the minimum thus is the potential evaluated at the classical mini- 
mum corrected by some terms, related to one-particle reducible diagrams and calculable 
in terms of potentials of lower order. The result for the 3-loop one-particle reducible 
contribution is 



1 



(V^2 + 



3 /, ^ 4 \ f2m 

2r+/?rw 



1 1 

log— + o 
mT 2 



+l3d\sc{h)mT 



9t 



+ 
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128 V 2 / (47r)3m2 



+ 



(A.16) 



where 



/?disc(^) 



9 

32 V 



(h + — \ 



-/i^-21)log3 + /i^log/i + 



+ (4 - 3/^2 + h^) log(2 + h)~ h\l + h^) log(l + h) + 

2 



2 /i2 4 



16 



2 + h. 



(A.17) 



Defining further 



/2(A) = y + A' 



2 ' 



(A.18) 
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2m 



9i 



(A.19) 



q2{h) = (2-/i2 + -/i4)log(2 + /i) 



2h 



+l/.Mog(3/.)-^log3-^-A-l/^^ 



(A.20) 



the vacuum energy density is 

m|(/x) 



V(v) = 



A, 



12% 



m|(6 + h^) + 



, t-2 - 2 "^93 



fx 



V my 2 



+ 



h A- — ^ /2m 

/3disc(/i)+/3(/i) 



rriT 2 



+ 52 
1 



+ 



+ 



9i 



(47r)= 



(A.21) 



Computing dV{v)/dm1 gives the final result for the scalar condensate in loop ex- 
pansion: 

3 



9l 



9l 



{k.22) 



2A3 



47r 
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/^V gl 



2A3 8(47r) 
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2^3 
9l 



f2{h)\og 



/i 



mT(/x) 



f'. 



2m 



9l 



6 



+ q2{h) 



rriT 



where the difi^erent loop contributions are on difi^erent lines and where (3{h) is the 3- 
loop quantity to be determined. Note how the general n dependence arises from an 
incomplete cancellation between the tree and 2-loop terms and how the /j, dependence 
of the 1-loop term is compensated to order by the 3-loop term. 



Appendix B 

In this appendix we describe an "improved" version of the 2-loop effective potential in 
Landau gauge used for the computation of the scalar condensate at different tempera- 
tures by the "Coleman- Weinberg" method. 
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It was pointed out in (see also Appendix A) that the effective potential in the 
Landau gauge contains singularities at the classical tree minimum, where the Goldstone 
mass 



ml = ml + A30^ (B.l^ 



vanishes. The order of leading singularities on the n-loop level can be found from power 
counting: 

sing ^ ^ (B.2) 



2 



where the C„ are some coefficients. These terms are not dangerous if the location of 
the minimum of the effective potential is far enough from the tree-level value (this 
is the Coleman- Weinberg regime). However, deep in the broken phase the quantum 
corrections are small, the Goldstone mass 1712 is close to zero and terms like ( |B.2| ) may 
become important. As we discussed in Sec. 5 of Q], due to the singular terms the 
computation of the ground state energy by the CW method, needed for the estimate 
of the value of the condensate (^V)) differs from the straightforward Ti expansion 
described in Appendix A by fractional powers of h\ at n loops, the difference is of 
order ^'^"'''2^ In particular, since the effective potential is known up to two loops, the 
difference between the two methods is of order . This is unacceptable since the 
^^-contribution to the effective potential is to be determined. 

Below we show how this disadvantage of the CW method can be removed, so that 
the CW method can be used both near and far from the critical temperature. 

The main idea is to redefine the 2-loop potential, including in it all leading singular- 
ities: 

00 

^Tor" = + E Vn'"''- (B.3) 

?i=3 

In this way, the expression for the ground state energy becomes analytic in as it 

must be in the broken phase, where there are no massless physical excitations. It is 

5 

interesting that the requirement that the ^2 correction be absent in the expression for 
the ground state energy is powerful enough to determine the structure of the leading 
singularities. 

Let us rewrite the tree and 1-loop effective potentials in the following form: 

\/o = ^/, V, = -^[Qml + m], + ?,z']-^mlz^ + 0{z'). (B.4) 



Here m^y = ^dsy {—'ml) / ^3 and mj^ = — 2m| are the tree-level values of the W and 
Higgs masses, z = m2 and 

ml='^glim^+^4mH). (B.5) 
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Then the solution to the minimization equation = has the form 



dz 

3 , 



z"^ = ml + 0{h^). (B.6) 

Now, if this solution is used in computing the values of the singular parts of the effective 
potential in eq. (|B.2|) , one can see that any of the n-loop contributions is of order . 
Since such a term must be absent, we get 

oo 

Y: Vr^iz' = ml) = 0. (B.7) 

n=l 

Therefore, the whole sum must have the form 

oo * 

YVr^{z) = --iz'-ml)l (B.8) 
in order to reproduce correctly the known 1-loop term ~ z^, which is non-analytic 

2 

m z . 

In this derivation the explicit form of the 2-loop potential was not used, so that 
the prediction of the term non-analytic in h"^ may be compared with the direct 2-loop 
computation. From 0, the non-analytic terms on the 2-loop level are 

indeed coinciding with the h expansion of eq. ( [B.8| ). 

To conclude, the improved form of the effective potential, which reproduces the 
correct structure of the h expansion of the ground state energy, has the form 

^^improved ^ y2-ioop(m2 ^ ml - ml) - —mlm^, (B.IO) 

where the last term must be subtracted in order to avoid double-counting. It is this 
potential which was used for the computation of the scalar condensate at fixed tem- 
perature by the CW method in Sec. H. 



Appendix C 

Here we calculate the W and Higgs correlators in the 3d theory of eq. ( p..l|) at 1-loop 
order. The Lagrangian masses of the vector, Higgs and Goldstone fields are in eq. ( |A.3|) . 
With these masses, the radiatively corrected propagators of the vector and Higgs fields 
are of the form 



k^A-v)A){v)) = 5^' ^, +;,/:''^^(p2^ + longitudinal part. 



-P)Mp)) = 2 TlH( 2V 
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To calculate the self-energies 11 , 11 , one needs the basic integrals 

1 m 



, 2 

pZ _|_ 

dp 



47r' 



[p2 + mj] [{p + ky + ml] 



mi + 1712 — i{k'^y^'^ 
7r(P)V2 "^mi + m2 + 2(A;2)i/2' 



In 



The integration measure here is 



d'^p 



(C.2) 



(C.3) 



(C.4) 



where d = 3 — 2e. 

The contributions of the diagrams in fig. |2^.a to the vector self-energy 11^ in eq. ( |U.1| ) 
are (V is a vector, S a scalar, and r] a ghost propagator; /c^ is the Euclidian external 
momentum) 



n 



w 

ss 



+Bo{k'^; ml, ml){k'^ + Amj 



k'^ + 2{ml + ml) + 



{m\ — miy-\ 



fc2 



+ 



2 2 

'mf — mo 



Ao(m?) - AoimD] (^^\^) - Ao(m?) - 3A,{ml) 



(C.5) 



n 



W 

BY 



n 



gl\Bo{P;ml,mT] 



5fc^ - 4m|, + 2 



+5o(A;2;m|,0) 



Am'^k'^ 
k^ 



k^ 

mJj. 



+5o(A;';0,0)(--^)+Ao(m|) 



12 ~ 2^ ^ 4A;2m| 



2 \2. 



-A;^ — 2m\ + 6m^ 



||so(A;';m2,m2^r .^2 , .^2 K ^t) 

+Eo(A;2;m2,0) 



A;2 



2\2 



+ ml) 



1 + 



fc2 

2 2 

m| — 



A;2 



2 

2\ 



F J' 



S'S l2 



, (C.6) 



(C.7) 
(C.8) 
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ttW 



9l 



(C.9) 



n 



(C.IO) 



For the 1-loop contributions to the Higgs self-energy 11^ in eq. (|C.1|) , one gets from 
the diagrams in fig. I23.b the results 



(C.ll) 



n 



H 
VV 



Bo{k ]mj,,mrj,) k + Ak rrij^ + Srurp 



-2Eo(fc2;m^,0)(A;2 + m^)2 
+Bo{k'^; 0, 0)k^ - 2Ao(m|)m| 



2 2 2 
T; '"2 



fc^ + 2P{m^ + ml) + (m^ - mj)'^ 



2\2 



+Ao(m; 



m^- — m2 



(C.12) 



(C.13) 



H 



n 



-3A, 



(C.14) 
(C.15) 
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